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THE  INSTITUTES  FOR  ENVIRONMENTAL  RESEARCH 


The  mission  of  the  Institutes  is  to  study  the  oceans,  and  in¬ 
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*  ^he  Institute  for  Oceanography 

works  to  increase  knowledge  and  improve 
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action  with  the  total  physical  environment 
of  the  globe. 

•  The  Institute  for  Atmospheric  Sciences 

seeks  the  understanding  of  atmospheric 
processes  and  phenomena  that  is  required 
to  improve  weather  forecasts  and  related 
services  and  to  modify  and  control  the 
weather. 
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Abstract 


y 

Methods  for  predicting  the  performance  of  high  frequency 
ionospheric  telecommunication  systems  are  revised  to  include  recent 
improvements  in  the  basic  ionospheric  and  geophysical  data.  Revisec 
techniques  for  the  processing  of  this  data  by  high  speed  electronic 
computers  is  described  in  detail  with  emphasis  on  better  statistical 
descriptions  of  the  expected  performance  of  radio  systems  depending 
upon  ionospheric  propagation  of  radio  waves.  The  application  of  the 
prediction  techniques  to  communication  problems  is  illustrated  and 
the  concept  of  service  probability  in  ionospheric  telecommunication 
systems  is  introduced.  Comparisons  between  predictions  and  circuit 
operations  are  shown. 
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1.  Introduction 

Recent  improvements  in  basic  ionospheric  and  geophysical  data  and 
the  increased  availability  of  electronic  computers  to  process  these  data 
have  made  it  possible  to  improve  previous  methods  of  predicting  the 
performance  of  ionospheric  telecommunication  circuits.  It  is  the  purpose 
of  this  report  to  outline  a  revised  prediction  method  based  on  an  assumed 
parabolic  distribution  of  the  electron  density  in  the  ionosphere. 

Much  has  been  written  on  parabolic  approximations  to  the  ionospheric 
layers,  ionospheric  radio  wave  propagation,  and  computer  programming. 
Specific  reference  is  made  to  a  number  of  articles  and  books  and  others 
are  listed  as  additional  references. 

It  is  not  the  intent  of  this  report  to  review  the  above,  but  to 
outline  a  method  of  predicting  the  performance  of  high  frequency  iono¬ 
spheric  telecommunications  systems  employing  those  methods  from  the 
literature  which  are  consistent  with  available  world-wide  data  and  which 
can  be  used  economically  while  maintaining  flexibility  to  easily  incor¬ 
porate  advances  in  ionospheric  research  and  data  collection  techniques. 


The  sections  to  follow  will  outline  (1)  philosophy  of  the  method, 

(2)  basic  ionospheric  data  available,  and  (3)  methods  for  predicting 
circuit  performance. 

2.  Philosophy  of  Prediction  Method 

The  need  for  high  frequency  prediction  methods  in  the  past  few  years 
has  resulted  in  the  development  of  many  diverse  models  to  represent  the 
factors  affecting  the  propagation  of  high  frequency  signals. 

The  models  range  from  very  simple  ones,  using  only  a  few  variables, 
to  very  elaborate  ray-tracing  techniques  which  require  a  precise  detailed 
knowledge  of  geophysical  and  ionospheric  parameters  to  yield  a  satisfac¬ 
tory  result. 

The  model  described  in  this  report  is  intended  to  be  as  sophis¬ 
ticated  as  the  basic  data  will  permit  and  is  designed  to  use  all  data 
available  on  a  world-wide  basis  to  predict  am  average  profile  of  electron 
density  versus  virtual  height  for  the  path  being  considered.  This  model 
retains  the  equivalence  theorem  and  transmission  curve  concept  to  remain 
consistent  with  the  methods  used  to  scale  and  predict  ionospheric  infor¬ 
mation.  The  electron  density  profile  along  the  path  is  assumedto  be 
adequately  represented  by  two  parabolic  layers.  The  height  of  maximum 
ionization.,  thickness  and  electron  density  are  derived  from  locations  near 
the  points  of  actual  reflection  along  the  path  instead  of  the  classical 
"two-control-point"  method  previously  used  in  the  calculation  of  the  upper 
limit  of  frequencies  and  transmission  loss. 

Averaging  of  geophysical  and  ionospheric  variables  is  effected  along 

the  path  to  yield  monthly  median  values  of  critical  frequencies  andlosses 

at  specific  frequencies.  These  values  are  combined  with  the  day-to-day 

distributions  to  predict  the  signal  level  exceeded  any  fraction  of  the 

days  within  the  month  (or  the  likelihood  of  a  given  level  being  exceeded 

on  a  randomlv  chosen  davl. 

/  /  *  • 
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In  this  report  the  convention  is  adopted  of  using  lower  case  letters 
to  denote  power  in  watts,  or  power  ratios  andcapital  letters  are  usedto 
denote  their  decibel  equivalents.  Thus  w  watts  may  be  expressed  as 
W  =  10  log  w  decibels  relative  to  one  watt.  This  convention  sometimes 

LO 

makes  it  impracticable  to  follow  the  recommendation  that  abbreviations 
for  units  based  directly  on  proper  names  be  capitalized.  However,  when 
no  confusion  will  arise  the  recommendation  will  be  followed,  e.  g.,  MHz 
(the  H  is  obviously  not  in  decibel  units). 

3.  Basic  Ionospheric  Data 

The  ionosphere  has  four  principal  regions,  or  layers,  which  affect 

9 

the  propagation  of  high  frequency  radio  waves.  They  are  the  D,  E,  FI, 
and  F2  layers  in  order  of  increasing  height  above  the  earth' s  surface  and 
increasing  ionization.  Either  the  E  or  F2  region  totally  reflect  radio  waves 
having  a  sufficiently  low  frequency  or  angle  of  incidence  while  most  of 
their  attenuation  occurs  during  their  passage  through  the  upper  D  region 
and  the  lower  E  region.  Basic  atmospheric  density  and  ionization  data  for 
these  main  layers  is  thus  required  and  used  in  this  model  to  predict  the 
strength  of  high  frequency  sky-wave  signals. 

3.  1.  F2 -Region  Ionization 

The  ionization  density  of  the  F2  region  is  very  unstable  and  exhibits 
marked  variations  with  geographic  latitude  and  longitude,  local  time, 
season  of  the  year,  and  solar  activity. 

The  Institute  for  Telecommunication  Sciences  aud  Aeronomy  predicts 
the  characteristics  of  the  F2  region  in  terms  of  the  vertical  incidence 
critical  frequency  (foF2)  and  the  M(3000)F2  factor  as  determined  from 
vertical  ionosonde  records  [Ostrow  1962], 

ThefoF2  may  be  scaled  directly  from  the  vertical  ionogram  (figure 
3.2).  Through  the  use  of  transmission  curves  [Smith  1939]  the  vertical 
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TIME  variation: 


ionogram  also  yields  the  M(3000)F2  which  relates  the  vertical  incident 
(foF2)  frequency  to  the  oblique  3000-km  frequency  as 

F2(3000)  MUF  =  foF2  •  M{3000)F2  (3.  D 

Tables  3,  1  and  3.  2  are  examples  of  monthly  median  world-wide 
numerical  coefficient  representation*  of  these  two  parameters  for  a  given 
month.  The  "numerical"  maps  represented  by  these  coefficients  denote  a 
function,  P{X,  8,  t),  of  the  three  varia.bles:  latitude  (X),  longitude  (9),  and 
local  time  (t).  The  function  P(X,  8,  t)  is  obtained  by  fitting  appropriate 
mathematical  functions  to  the  observed  ionospheric  sounding  data  collected 
from  a  world-wide  net  of  vertical  sounders. 

The  general  form  of  P(X,  9,  t)  is  the  Fourier  time  series 

H 

r(X,  9,  t)  =  a0  (X,  9)  +  Z  [a^  (X,  8)  cos  j  t  +  b^  (X,  0)  sin  j  t]  (3.2) 

j  =  1 

LOstrow  1962]. 

H  denotes  the  number  of  harmonics  retained  to  represent  the  diurnal 
variation.  The  Fourier  coefficients,  a^  (X,  0)  and  b5  (X,  6),  which  vary 
with  the  geographic  location,  are  represented  by  a  series 

K 

£  D„k  Gk  (X,  9),  (3.3) 

k=0 

where  the  Gk  (X,  0)  are  shown  in  table  3.3.  The  index  (s)  denotes  which 
Fourier  coefficient  is  represented,  in  the  order 

D<>k  =  Fourier  coefficients  defining  F(X,  8,  t) 


s  =  2  j,  for  aj  (X,  0),  j  =  0,  1,  2 . H 

s  =  2  j  -  1  for  b3  (X,  9),  j  =  1,  2,  3 . H 


The  above  function  (T)  has  the  same  form  for  both  the  foFJ  and  the 
M{3000)F2  which  are  utilized  by  this  method  [Jones  and  Gallet  1962a].  A 

revised  method  for  the  generation  of  foF2  and  M(3000)F2  will  be  used  in 
the  near  future  [Jones,  et  al  1966]. 
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3.2.  FI -Region  Ionization 

The  FI  layer  also  affects  HF  propagation  especially  in  the  daytime. 
However  the  geographic,  time  and  solar  activity  dependence  described 
above  for  the  F2-layer  parameters  is  not  presently  available  for  the 
FI  layer.  This  layer  is  therefore  not  included  as  a  separate  layer  in  the 
present  method  although  the  effect  of  FI -layer  retardation  on  waves 
propagated  via  the  F2  layer  is  partially  accounted  for  (see  section  3.4.  on 
layer  heights), 

3.3.  Regular  E-Region  Ionization 

The  regular  E  layer  does  not  exhibit  irregularities  as  complex  as 
those  associated  with  the  F2  region  since  its  ionization  is  controlled 
primarily  by  the  sun' s  zenith  angle.  Based  on  the  Chapman  theory  and 
using  the  best  current  estimates  of  the  rates  of  electron  production  and 
loss,  the  following  semiempirical  equation  has  been  derived  (Knecht  1962), 

1 

foE  =  0.9  [(180  +  1.44  5)  cos  HU4  (3.4) 

S  =  12  month  running  Zurich  sunspot  number 

where  r  , 

Hi  =  zenith  angle  of  the  sun  -  degrees 

foE  =  monthly  median  critical  frequency  of  E  layer  -  MHz. 

Although  the  above  equation  adequately  lepresents  the  E-region 
electron  density  during  most  daylight  hours,  it  fails  to  predict  E-layer 
critical  frequencies  during  the  twilight  hours  (H<>  90°)  since  an  foE  of 
zero  at  a  sun' s  zenith  angle  of  90  degrees  is  clearly  unreasonable,  e.g.  , 
E-layer  critical  frequencies  below  700  KHz  have  seldom  been  observed 
(Watts  and  Brown  1954). 

Since  the  D-region  ion’ospheric  absorption  is  closely  correlated 
with  the  '-ritical  frequency  of  the  E  region,  it  is  convenient  to  relate  foE 
to  the  index  of  ionospheric  absorption  which  is  also  a  function  of  solar 
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activity  and  the  zenith  angle  of  the  sun.  The  ionospheric  absorption  index 
is  expressed  as 

I  =  (1  +  .  0037  S)  cos  (K  (3.  5) 

S  =  12  month  running  average  Zurich  sunspot  number 

where 

=  sun1  s  zenith  angle  -  degrees 

K  =  90°/  102°  =  0.  881  (102°  is  sunset  at  110  km 
height  of  the  E  layer). 

The  above  value  of  I  is  related  to  the  E(2000)MUF  by  figure  3.  1 
which  is  mathematically  represented  by  a  polynomial: 

E(2000)  MUF  =  3. 41  +  38.  43  •  I  -  68.  07  •  I2  +  89.  97  •  I3  (3.  6) 
-70.  97  •  I4  +  29.  51  •  Is  -  4.  99  •  I6  (MHz) 

The  above  equation  yields  an  E(2000)MUF  of  3.41  MHz  during  the 
nighttime  hours  corresponding  to  an  foE  of  0,  7  MHz  (Watts  and  Brown 
1954).  The  regular  E  layer  is  assumed  to  be  very  predictable  and  its 
associated  distribution  is  considered  negligible. 

3.  4„  Height  of  Maximum  Ionization  of  the  F2  and  E-  Layers 

The  height  of  maximum  ionization  of  the  F2  region  (hmF2)  is 
needed  in  practical  communication  problems  and  geophysical  research. 
Unfortunately,  direct  measurement  of  hm  is  impossible  from  vertical 
incidence  ionograms.  Although  accurate  computer  methods  of  determining 
electron  density  profiles  from  virtual  height  curves  are  available  [Wright 
and  Paul  1963],  this  information  is  not  currently  available  on  a  world-wide 
basis  and  another  way  of  estimating  the  hm  of  the  F2  region  is  required. 

As  an  interim  method  it  is  proposed  to  use  the  existing  world-wide 
numerical  maps  of  M(3000)F2  and  a  linear  conversion  formula  developed 
by  Shimazaki  [  1 9 S 5 .]  to  obtain  the  necessary  estimates  of  hm.  There  has 
been  some  criticism  of  Shimazaki1  s  formula  [Wright  and  McDuffie  I960] 
but  its  use  in  practical  prediction  methods  is  considered  to  be  justified 
as  is  explained  later  in  the  discussion. 
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ABSORPTION  INDEX 


In  practical  prediction  methods,  such  as  this,  it  is  not  always 
possible  to  justify  rigorously  every  technique  and  parameter  employed, 
e.g.  ,  it  should  be  restated  that  transmission  via  the  FI  region  on  oblique 
paths  is  not  considered  in  this  method;  however,  some  adjustment  to 
h^FZ  needs  to  be  considered  when  using  vertically  incident  ionograms 
which  include  FI  layer  retardation.  Improvements  depend  in  a  large  part 
on  the  improvements  in  basic  data,  scaling  methods,  and  basic  research. 
Additional  background  for  this  material  can  be  found  in  Davies  [1965]. 

The  maximum  usable  frequency  factor,  M(3000)F2,  for  propagation 
by  way  of  the  F2  layer  over  a  3000-kilometer  path  is  regularly  scaled 
from  ionograms  on  a  world-wide  basis.  Another  parameter  sometimes 
scaled  regularly  is  hpF2,  the  virtual  height  on  the  ionogram  at  a 
frequency  equal  to  0.  834  foF2.  This  frequency  was  selected  for  scaling 
because  in  parabolic  layer  theory,  the  virtual  height  at  83.4  %  of 
the  critical  frequency  of  the  layer  is  equal  to  the  true  height  (h^)  of 
maximum  ionization  (see  figure  3.2).  Shimazaki  [  1955]  has  demonstrated 
that  these  two  parameters  are  approximately  related  by  the  foliowing 
linear  equation: 

h  F2  3s  -  176  4  1490  (3.7) 

P  M(3000)F2 

This  formula  was  derived  theoretically  using  plausible  approxi¬ 
mations  for  the  layer  shape.  Shimazaki  also  compared  values  of  h  F2 

lr 

calculated  from  M(3000)F2  by  the  formula  with  the  corresponding  scaled 
values  of  hpF2  for  18  stations  for  the  months  of  March,  June,  September, 
and  December  1932,  and  for  Tokyo  for  the  same  months  and  the  years 
1949  through  1954  (see  figure  3.3).  The  daily  mean  difference  between  the 
two  values  rarely  exceeded  6.  0  %  ,  thus  showing  good  agreement  for  all 
seasons,  geographic  locations  and  levels  of  solar  activity. 
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INCREASING  FREGUENCV 


The  MUF  factor,  M(3000)F2,  is  scaled  from  the  ionograms  along 
with  foF2,  and  both  are  mapped  using  a  numerical  method  on  a  world-wide 
basis.  Since  hpF2  is  not  mapped  world-wide,  the  present  prediction 
program  finds  it  from  Shimazaki' s  formula  (3.  7). 

In  I960  Wright  and  McDuffie  reconsidered  Shimazaki' s  formula 
(along  with  other  methods  of  inferring  the  height  of  the  F-region  peak) 
in  the  light  of  h^  values  obtained  from  the  computer  methods  of  deriving 
N(h)  profiles  from  ionograms.  Their  analysis  showed  that  h  F2  was  "free 
from  bias  at  night  at  low  and  medium  latitudes,  but  that  systematic  differ¬ 
ences  are  significant  in  the  daytime  and  at  high  latitudes  at  all  times.  " 

For  these  latter  times  and  latitudes,  h^  is  significantly  less  than  hpF2 
as  deduced  from  Shimazaki1  s  formula,  but  an  empirical  linear  relationship 
between  h^  and  M(3000)F2  which  fits  the  data  quite  well  can  usually  be 
found  though  the  slope  and  intercept  values  may  be  a  function  of  time, 
location,  etc. 

Wright  [1964]  later  indicated  that  hpF2  as  inferred  from  M(3000)F2 
would  be  preferable  to  for  prediction  purposes  even  at  che  times  when 
hpF2  is  greater.  This  is  because  the  M(3000)F2"»hp  theory  does  not  take 
into  account  the  effects  of  the  daytime  E  and  FI  layers  on  M(3000)F2.  The 
virtual  height  of  the  observed  F2  trace  includes  the  effects  of  retardation 
in  the  lower  layers.  Therefore  the  application  of  the  3000-km  transmission 
curve  to  such  an  ionogram  results  in  an  M(3000)F2  smaller  than  that  for  the 
F2  layer  alone.  The  smaller  M(3000)F2  yields  a  greater  hpF2  than  would 
be  obtained  for  the  F2  layer  alone.  While  this  effect  is  a  disadvantage  for 
deducing  true  heights  for  geophysical  research,  it  is  an  advantage  in 
computing  elevation  angles  for  communication  paths  via  the  ^2  layer  when 
independent  information  on  the  FI  layer  parameters  is  not  otherwise 
available. 

Figure  3.2  illustrates  how  changes  in  the  virtual  height  (vertical 
scale)  and  semithickness  of  the  F2  layer  affect  the  values  of  M{3000)F/. 
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factors  scaled  from  ionograms.  The  thick  solid  line  is  the  3000-kilometer 
transmission  curve  with  MUF  lector  marked  on  a  logarithmic  horizontal 
scale.  The  MUF  factor  increases  with  decreasing  values  of  the  virtual 
height,  but  must  by  definition  be  read  where  the  h^f  curve  becomes 
vertical,  i.  e.  ,  at  the  critical  frequency  (long  vertical  lines).  Note  that 
the  factor  is  not  read  at  the  point  of  tangency  to  the  transmission  curve. 

The  short  vertical  lines  on  the  h/~f  curves  correspond  to  83.  4%  of  the 
critical  frequency  and  thus  intersect  the  h.'  -f  curves  at  a  virtual  height 
equal  to  the  true  height  of  maximum,  1^,  of  a  single  parabolic  l^yer. 

Note  that  hj^  is  not  necessarily  equal  to  the  height  of  tangency  of  the 
transmission  curve,  h^an  (crosses). 

Since  the  frequency  scale  of  the  ionogram  is  logarithmic,  the  MUF 
factor  scaled  is  independent  of  the  absolute  value  of  the  critical  frequency 
(at  least  as  fa.r  as  the  scaling  procedure  is  concerned).  The  two  thin  solid 
lines  are  h^f  curves  for  parabolic  layers  of  equal  semithickness  (100  km) 
but  at  different  heights  illustrating  the  increase  of  MUF  factor  with  decreas¬ 
ing  height.  The  increase  of  MUF  factor  with  decreasing  semithickness  is 
illustrated  by  the  dotted  h7-f  curve  which  is  tangent  at  the  same  virtual 
height  (crosses)  but  has  a  semithickness  of  35  Ion  instead  of  100  km. 

Since  scaling  the  MUF  factor  is  strictly  a  mechanical  procedure, 
it  makes  no  distinction  between  increases  in  the  true  height  of  the  F2  layer 
and  apparent  height  increases  due  to  the  retardation  in  the  FI  layer  beneath. 
Therefore,  when  the  FI  layer  is  present  during  the  daytime  hours,  neither 
MUF  factors  nor  values  of  hpF2  scaled  from  the  same  ionogram  (at  0.  834 
foF?,)  will  correspond  to  the  true  height  of  the  F2  layer.  However,  they 
may  still  be  linearly  related  to  a  good  approximation  through  Shimazaki' s 
formulas  as  his  figures  demonstrate. 

Figure  3.4  illustrates  the  relationship  between  an  observed  ionogram 
(solid  line)  when  both  FI  and  F2  layers  are  present  and  an  ionogram 
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calculated  from  the  true  height,  h^FZ,  of  the  F2  layer  (dotted  curved 

line).  Portions  of  the  3000-kilometer  transmission  curves  (marked  A  and 

B)  are  shown  tangent  to  each  ionogram.  The  M(3000)F2  factors  may  be 

read  from  the  scales  at  the  critical  frequency  of  the  F2  layer.  Note  that 

the  MUF  factor  (scale  A)  for  the  observed  ionogram  is  less  than  that 

(scale  B)  for  the  h/-f  curve  corresponding  to  h  F2  (marked  h/F2  only). 

m 

In  this  case  hi  F2  was  350  kilometers,  while  the  values  of  h  F2 
m  p 

obtained  by  Shimazaki' s  formula  were  341  km  for  the  larger  MUF  factor 
(scale  B)  and  360  km  for  the  small  (scale  A).  Thus  a  h  -f  curve  corres¬ 
ponding  to  h^_F2  =  360  would  everywhere  be  10  km  higher  than  the  dotted 
curve  and  would, therefore, fit  the  observed  ionogram  somewhat  better  at 
least  for  frequencies  above  the  critical  frequency  of  the  FI  layer. 

Virtual  heights  found  in  the  above  manner  should,  therefore,  yield 
more  accurate  elevation  angles  (by  means  of  the  iterative  process  discussed 
elsewhere  in  this  report)  than  would  virtual  heights  based  on  the  actual 
true  heights  from  which  the  FI  retardation  had  been  removed.  Of  course, 
if  the  FI  layer  were  explicitly  included  as  a  separate  layer,  the  unretarded 
true  F2-layer  heights  would  be  preferred. 

The  E  region  is  considered  as  a  separate  layer  in  this  method;  thus 
some  allowance  is  made  for  the  effect  of  the  retardation  in  the  E  region 

upon  the  calculated  value  of  h  F2  derived  from  the  M(3000)F2  scaled  at 

P 

vertical  incidence.  The  formula  to  follow  is  a  general  formula  which  could 
be  used  to  calculate  the  retardation  in  any  underlying  parabolic  layer. 

The  amount  of  retardation  expected  when  a  vertical  ray  penetrates 
the  E  region  (f  >  fc )  assuming  a  parabolic  distribution  of  electron  density 
with  height  is 

-Ah  =  {£/£.  log[f^44]-2}-y 

f  =  frequency  at  which  hpF2  was  scaled  (=  0.  834  foF2) 
where  fc  =  critical  frequency  of  the  E  region  (feE) 

yE  =  semithickness  of  E  region. 
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0.834foF2/foE 

Figure  3.  5.  Correction  to  h  F2  Due  to  Retardation 
in  E -Region  (Vertical  Incidence) 
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Figure  3.  5  shows  some  typical  values  of  vertical  retardation  as  a 
function  of  foF2/  fo(E).  A  complete  derivation  of  the  above  relationship 
is  available  [Kelso  1964],  in  terms  of  group  path  and  the  parabolic  layer 
assumption. 

The  maximum  height  of  the  E  region,  hm(E),  is  taken  to  be  constant 
for  all  hours  and  seasons  at  130  km  [Al'pert  I960], 

3.  5.  Quasi-Empirical  Determination  of  F2-and  E-Layer  Semithickness  (y) 

The  height  of  maximum  ionization  of  the  F2  layer  hmF2  in  units  of 
the  semithickness  y  is  mapped  as  a  function  of  geomagnetic  latitude  and 
solar  zenith  angle  (figure  3.  6  and  3.  7)  for  high  and  low  solar  activity.. 

The  contour  plots  were  generated  by  using  smoothed  values  of  the  height 
of  maximum  electron  density  and  quarter -thickness  (SCAT  =  y /  2)  from 
observations  primarily  along  the  75th  meridian  [Wright,  Wescott,  and 
Brown  I960]. 

June  and  January  data  were  used  to  estimate  extreme  values  of  the 
zenith  angle  of  the  sun  at  specific  geographic  latitudes.  A  scatter  plot 
of  1^/ SCAT  versus  sun' s  zenith  angle  produced  average  regression  lines 
with  little  scatter. 

Stations  off  the  75th  meridian,  varying  in  geomagnetic  latitude, 
were  checked  against  values  on  the  75th  meridian  to  verify  the  relationship 
with  geomagnetic  latitude.  The  values  compared  favorably,  thus  geo¬ 
graphic  longitudes  and  latitudes  were  replaced  by  geomagnetic  latitudes. 
The  ratio  (^m/y)  depends  somewhat  on  solar  activity.  Therefore  contour 
charts  were  prepared  for  high  (125)  and  low  (25)  twelve  month  running 
average  Zurich  sunspot  numbers  [Lucas  and  Remmler  1966]. 

The  ratio  h^/  y  is  overestimated  when  hpF2  is  used  as  an  approx¬ 
imation  to  hmF2,  especially  during  the  daylight  hours  in  the  summertime. 
The  reasons  for  this  are  as  explained  in  the  preceding  section.  However, 
until  the  FI  layer  is  represented  separately,  the  values  of  semithickness, 
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Figure  3.  7.  Contour  Chart  of  the  Ratio  hmF2/  yF2 
for  Low  Solar  Activity  (SSN  =  2  5) 


using  hpF2  as  an  estimate  of  l^FE.  yield  a  more  realistic  description 
of  the  complete  F  layer  than  values  estimated  from  the  true  height  of  the 
F2  region.  This  is  especially  true  in  a  practical  prediction  routine  that 
relies  upon  the  equivalence  theorem  and  the  sec  0  "corrected"  law. 

The  semithickness  of  the  E  region  is  taken  to  be  stable  throughout 
the  day  and  for  all  seasons  as  20  km  [Al'pert  I960]. 

4.  Basic  Noise  Data 

4.  1.  Median  and  Decile  Values  of  Man-Made  Noise 
The  available  man-made  noise  power  at  the  input  of  an  equivalent 
loss-free  receiving  antenna  used  in  this  model  is  assumed  to  vary  as 

N,  =  V  +  12.60  log,  (f/3) 

Nb  =  man-made  noise  power  in  1  Hz  band  - 

,  dB  <  1  watt 

where 

V  =  measured  value  of  man-made  noise  power 
at  3  MHz  in  1  Hz  band-dB  <1  watt 

f  -  operating  frequency  (3-30  MHz) 

If  measured  values  of  man-made  noise  are  not  available,  an  estimate 
may  be  obtained  from  the  population  of  the  area  about  the  receiving  site 
as  shown  in  figure  4.  1.  The  upper  and  lower  decile  values,  Du  and  D^, 
are  taken  to  be  9  dB  and  7  dB  from  the  graphical  values  which  are  the 
medians,  respectively  [Spaulding  1965], 

4.2.  Median  and  Decile  Value  of  Galactic  Noise 
The  galactic  noise  estimates  and  the  associated  distributions  are 
those  levels  extrapolaced  to  3  MHz  from  Cottony  and  Johler  [1952],  and 
verified  using  a  vertical  antenna  [Crichlow  and  Spaulding  1965], 

The  median  galactic  noise  is  represented  by 
Ng  =  165  +  9.  555  log,  (f/3) 

where  Nt  -  expected  median  value  of  the  galactic  noise  in  a 

1  Hz  band  -  dB  <  1  watt 

f  =  operating  frequency  -  MHz  . 


22 


1 


The  variability  (Du  and  D^)  is  taken  to  be  2  dB  about  the  median 
[CCIR  322]  at  times  when  it  is  possible  to  see  the  galaxy  at  the  operating 
frequency.  The  ionosphere  above  the  receiving  locations  is  examined  to 
determine  whether  ionospheric  penetration  is  likely  and  galactic  noise 
present. 

4.3.  Median  and  Decile  Value  of  Atmospheric  Radio  Noise  (3-30  MHz) 

The  values  of  atmospheric  radio  noise  used  in  these  methods  are 
taken  from  the  world-wide  1  MHz  noise  maps,  frequency  dependencies 
(3-30  MHz),  and  variability  charts  (3-30  MHz),  found  in  CCIR  Report  322. 

The  world-wide  (1  MHz)  atmospheric  noise  maps  are  for  discrete 
four-hour  time  blocks  for  four  seasons  of  the  year.  The  values  from 
two  adjacent  maps  are  used  and  an  interpolation  made  on  time.  A 
frequency  dependence  and  variability  accompanies  each  four -hour  time 
block  and  adjacent  values  are  used  to  interpolate  for  a  given  time. 

Methods  of  mapping  the  noise  and  the  accuracy  to  which  it  was 
mapped  are  explained  by  Lucas  and  Harper  [1965].  A  sample  numerical 
world-wide  map,  frequency  dependence,  and  variability  chart  appear  in 
figures  4.  2  and  4.  3  (k,  Boltzmann' s  constant;  t,  288°  Kelvin;  b,  effective 
noise  bandwidth). 

The  numerical  coefficients  of  table  4.  1  which  describe  the  world-wide 
distribution  of  the  1  MHz  atmospheric  noise  (dB  >  ktb)  are  evaluated  by 
the  function 

F  (X,  9)  =  {  Y,  [  C?  V*  sin  "  0  )  +  X  ]  sin  K  *  }  +  I  ;  (4.3) 

K 

where  X  =  normalizing  coefficient  in  longitude 

§  =  a  +  3  X, 
a  =  F  (0,  9) 

0  =  F  (tt,  0)  -  F  (0,  9) 

TT 

9  =  geographic  latitude  -  degrees  (north  positive,  south 
negative) 

X  =  geographic  longitude  -  degrees  (east  of  Greenwich). 
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The  coefficients  x,  a.  and  8  appear  on  the  top  of  tahie  4.  1.  The 
seen  in  the  above  table  are  in  ascending  order  in  k. 

The  "numerical"  representation  of  the  frequency  dependence  of  the 
atmospheric  noise  (figure  4.  3)  is  evaluated  by 

Y  (X,  N)  =  A-l  (N)  +  As  (N)  X  +  Aa  (N)  X2  +  .  .  .+  A7  (N)  X6  (4.  4) 
where 

At  (N)  =  b1(1  +  b1>2  N,  i  =  1,  .  .  .  ,  7 

N  =  amplitude  of  the  hourly  median  1  MHz  atmospheric 
radio  noise  (dB  >  ktb),  and 

X  =  8  x  2loglof  -11 
4 

where  f  is  operating  frequency  -  MHz. 

The  "numerical"  coefficients  shown  in  table  4.2  representing  tl 
day-to-day  distribution  (Du  &  D^)  of  the  atmospheric  radio  noise  are 
evaluated  by 

Y  (X)  =  A1  +  As  X  +  A3  X2  +  A4  X3  +  As  X4  ,  (4.  5) 

where 

X  =  log  (f),  and 

f  =  operating  frequency  -  MHz  . 

All  coefficients  shown  in  tables  4.  1  and  4.  2  are  listed  in  ascending 
order  from  left  to  right. 

4.  4.  Combination  of  the  Noise 

The  atmospheric,  man-made,  and  galactic  noises  are  evaluated  at 
the  receiving  site  and  the  predominant  median  noise  and  its  associated 
distribution  taken  to  be  the  prevailing  noise  power  at  all  probability 
levels.  Due  to  the  uncertainties  associated  with  the  predicted  values 
of  noise, a  negligible  additional  error  is  anticipated  through  the  use  of 
this  assumption. 
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Fourier  Coefficients  Representing  the  1  MHz  World-wide 
Distribution  of  Atmospheric  Radio  Noise, 
December-January-February  (0000-0400  Local  Mean  Time) 
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5.  General  Problem  of  Predicting 
Sky-Wave  Radio  System  Performance 

The  performance  of  ionospheric  telecommunication  circuits  varies 
markedly  with  frequency,  geographic  location,  time  of  operation,  equip¬ 
ment  parameters,  etc.  Much  of  this  variation  is  directly  related  to 
changing  ionospheric  conditions  which  (among  other  things)  dictate  the 
likelihood  of  long-distance  communications.  The  changes  affect  the 
attenuation  of  the  radio  signal  and  control  the  vertical  angles  of  arrival 
and  departure. 

Predictions  of  ionospheric  telecommunication  system  performance 
can  be  expressed  in  many  ways.  This  report  expresses  the  expected 
performance  as  follows: 

(1)  the  probability  that  a  sky-wave  path  exists  for  a  given 
frequency  (qf)  ; 

(2)  the  maximum  frequency  corresponding  to  a  specified 
probability  of  propagation  MUF  (qf) ; 

(3)  the  loss  in  the  communication  system  (La )  ; 

(4)  available  median  signal  at  the  receiving  system  input  (5); 

(5)  available  signal-to-noise  ratio  at  the  receiver  (s/n); 

(6)  likelihood  that  a  specified  signal-to-noise  ratio  will  be 
exceeded  (q^); 

(7)  probability  that  a  specified  grade  of  service  (gr)  will  be 
equalled  or  exceeded  for  a  specified  fraction  of  days  within 
the  month  (service  probability,  Or) [Bar sis  et  al  1961;  Norton 
1962;  Rice  et  al,  1965];  and 

(8)  the  minimum  operating  frequency  for  a  specified  fraction 
of  days  within  the  month  (LUF). 

5.  1.  Likelihood  of  a  Sky-Wave  Path 

The  detailed  considerations  in  the  determination  of  sky-wave  paths 
and  the  associated  probability  of  the  existence  of  a  sky-wave  path  are 
given  in  section  10.  In  general,  the  probabilitv  is  estimated  as  a  function 
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of  frequency,  for  a  specified  path  length,  at  a  given  geographic  location 
and  time.  The  probability  depends  upon  the  maximum  ionisation  level 
and  the  corresponding  height  in  the  various  ionospheric  regions.  Prob¬ 
ability  predictions  involve  both  estimates  of  median  ionospheric  conditions 
and  estimates  of  the  distribution  about  the  median  and  are  expressed  as 
the  fraction  of  days  within  the  month  that  some  sky-wave  path  is  expected 
to  exist  at  a  given  hour. 

5.2,  The  Maximum  Frequency  Corresponding  to  a 
Specified  Probability  of  Propagation 

The  maximum  frequency  corresponding  to  a  specified  probability 
of  propagation  involves  the  same  considerations  outlined  in  section  5.  1 
above,  except  that  the  probability  is  fixed  and  the  corresponding  frequency 
is  determined  rather  than  the  frequency  being  fixed  and  the  probability 
determined.  Normally  the  probability  is  fixed  at  0.  5  and  the  frequency 
corresponding  to  this  probability  is  referred  to  as  the  monthly  median 
Maximum  Usable  Frequency  or  MUF  (0.  5).  Note  that  MUF(q)  is  the 
quantile  of  order  q  of  the  distribution,  i.  e.,  MUF  <MUF(q)  with  prob¬ 
ability  q.  As  q  varies  from  0  to  1  the  quantile  of  order  q,  MUF(q),  will 
increase  from  zero  to  its  maximum  value. 

5.3.  The  Loss  in  the  Communication  System  (L3) 

System  loss  for  a  particular  propagation  path  in  high  frequency 
communication  circuits  may  be  defined  as  follows: 

L8  =  Lw  +  Lt  +  L*  -  Gt  -  Gr  +  Yp,  (dB)  (5.  1} 

where: 

Ls  =  system  loss  (signal  power  available  at  the  receiving 
antenna  terminals  relative  to  that  available  at  the 
transmitting  antenna  terminals  in  decibels) 

LbJ.  =  free  space  transmission  loss  based  on  ray  path  distance 

of  the  path  being  considered  and  the  radio  frequency  (decibels) 

Lj  =  losses  due  to  ionospheric  absorption  (decibels) 
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Ls  =  losses  due  to  ground  reflection  (decibels) 

Gt  and  Gr  =  the  transmitting  antenna  and  receiving  antenna  power 
gains  in  decibels  relative  to  an  isotropic  antenna  £ 

Yp  =  loss  as  so;  .ated  with  day-to-day  variations  in 

ionospheric  and  circuit  parameters  (statistically 
determined  and  dependent  upon  circuit  location, 
time,  path  length,  and  fraction  of  days  being 
considered). 

In  high  frequency  communication  circuits,  several  propagation 
paths  are  often  pos  !ible,  e.  g.  ,  a  single  reflection  from  the  F  region 
(IF);  a  single  reflection  from  the  E  region  (IE),  multiple  reflection 
from  the  E  and  F  region  (2 F,  3F,  2E,  etc.);  or  paths  involving  reflection 
from  both  regions  (1E1F,  1F1E,  etc.).  The  probable  paths  are  dependent 
upon  the  geometry  of  the  circuit  involving  layer  heights  and  great  circle 
distance  and  upon  the  relative  ionization  in  the  various  regions. 

For  most  systems  applications,  it  is  considered  adequate:  (1)  to 
evaluate  ,  L,  ,  ,  Gt  and  Gr  for  each  of  the  likely  propagation  paths, 

(2)  to  select  the  path  with  the  minimum  loss,  and  (3)  to  combine  this 
loss  with  the  empirically  determined  (Yp  )  which  includes  the  effect  of  the 
day-to-day  variations  in  the  parameters  used  in  estimating  Lbf,  , 

Gt  and  Gp  ,  plus  such  factors  as  ionospheric  focusing,  deviative  abso  •  'ion, 
polarization  mismatch  at  the  receiving  antenna,  and  the  addition  of  signals 
via  various  paths. 

5.4.  Available  Signal  at  the  Receiver 

The  available  signal  level  at  the  receiver  requires  only  that  the  loss 
in  the  communication  system  as  defined  in  section  5.  3.  be  combined  with 
the  power  available  at  the  transmitting  antenna  terminals.  The  fraction 
of  time  a  given  signal  level  is  exceeded  is  directly  related  to  the  fraction 
of  time  associated  with  Yp  in  section  5.  3. 

*  In.  this  report  Gt  and  Gf  are  in  the  direction  of  the  propagation  path 
and  include  all  losses  so  that  Gt  f  Gr  approximates  the  path  antenna 
goin  Gp  [ Rice,  et  al  1965]. 
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5.  5.  Available  Signal-to-Noise  Ratio  at  the  Receiver 
An  estimation  of  the  available  signal-to-noise  ratio  requires  that 
the  expected  available  signal  power  be  combined  with  the  expected  noise 
power.  The  expected  noise  power  includes  estimates  of  atmospheric, 
cosmic,  and  man-made  noise  powers.  Distributions  as  well  as  median 
values  of  both  signal  and  noise  are  involved  including  seasonal,  geograph¬ 
ic,  and  frequency  dependence.  Normally  the  expected  signal-to-noise 
ratio  is  expressed  as  a  median  of  the  hourly  medians  for  those  days 
some  sky-wave  path  exists  within  the  month. 

5.  6.  Likelihood  that  a  Specified  Signal-to-Noise  Ratio 

Will  be  Exceeded 

The  likelihood  that  a  specified  signal-to-noise  ratio  will  be  exceeded 
involves  the  same  considerations  as  in  section  5.3  except  that  signal-to- 
noise  ratio  is  fixed  and  the  likelihood  determined  plus  a  consideration  of 
the  probability  that  some  sky-wave  path  exists.  The  likelihood  that  a 
specified  signal-to-noise  ratio  is  available  at  the  receiver  input  is  referre 
to  as  circuit  reliability  q(Rh )  and  is  obtained  as  a  product  of  the  likelihood 
that  a  given  signal-to-noise  ratio  exists  on  the  days  a  sky-wave  path  is 
expected  (qs/n)  with  the  probability  that  a  sky-wave  path  exists  (qf ). 

5.  7.  The  Minimum  Operating  Frequency 
As  the  operating  frequency  : s  decreased,  the  likelihood  that  a 
specified  signal-to-noise  ratio  is  equalled  or  exceeded  also  decreases. 

The  minimum  operating  frequency,  or  Lowest  Useful  Frequency  (LUF), 
is  therefore  the  lowest  frequency  which  can  be  expected  to  have  a  signal- 
to-noise  ratio  equal  to  or  greater  than  that  required  for  at  least  a  given 
fraction  of  days  (qr )  within  the  month  (T).  The  required  signal-to-noise 
ratio  (Rh)  is  associated  with  the  type  and  quality  of  communication 
required  (grade  of  service,  gP)  and  is  expressed  as  an  hourly  median 
signal-to-noise  ratio.  Unless  otherwise  specified  the  LUF  is  estimated 
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on  the  basis  that  the  required  signal-to-noise  ratio  be  equalled  or  exceeded 
on  0.  90  of  the  days  of  the  month. 

6.  Theoretical  Basis  for  Communication  Path  Geometry 
6.  1.  Brief  Description  of  Model 
The  propagation  prediction  model  relies  primarily  upon  vertical 
incidence  ionosonde  data  to  predict  the  performance  on  oblique  incident 
communication  paths.  Probing  at  a  given  frequency  through  assumed 
parabolic  layers  (E  and  F2)  is  used  to  predict  the  virtual  height  of  reflec¬ 
tion  on  oblique  paths.  The  data  used  for  the  above  approximation  are  the 
monthly  median  foE,  foF2,  M-3000  factor,  F2  and  E  semithickness,  and 
probing  frequency  (f). 

6.2.  Vertical-to-Oblique  Transformation 
Since  the  basic  data  available  on  a  world-wide  basis  are  obtained 
from  vertical  ionosondes,  a  transformation  must  be  effected  in  the  appli¬ 
cation  to  an  oblique  path.  The  relation  between  the  frequency  (f)  of  an 
oblique  incident  wave  reflected  from  the  ionosphere  and  the  frequency {i') 
of  a  vertical  incident  wave  reflected  from  the  same  region  is  the  point 
of  interest.  Since  the  two  waves  are  reflected  from  the  same  virtual 
height,  the  relative  electron  density  at  the  level  of  reflection  must  be  the 
same  for  both. 

Disregarding  the  terms  due  to  the  earth' s  magnetic  field  and  also 
using  Snell' s  law  applied  to  a  flat  ionosphere  for  the  refractive  index  of 
the  ionosphere,  the  following  relation  is  obtained: 

f  =  f ,  sec  0,  (6.  1) 

where 

0  =  incident  angle  of  the  wave  upon  the  ionosphere 
f  =  oblique  frequency  at  h/ 
f vertical  incident  frequency  at  h/ 
h  =  virtual  height  of  reflection. 


For  the  curved  ionosphere  the  above  relationship  must  be  corrected  to 

f  =  f /  k  sec  0,  (6.  2) 

where  k  sec  0  is  referred  to  as  secant  0  (corrected).  The  ionospheric 
curvature  correction  factor,  k,  is  shown  as  a  function  of  the  transmission 
distance  in  figure  6,.  1  and  discussed  in  section  7.  1. 

Figure  6.  2  shows  the  vertical-to-oblique  transformation  graphically. 
The  dashed  curve  C,  is  the  so-called  transmission  curve  (given  distance) 
by  which  vertical  ionograms  are  scaled.  The  point  at  which  the  trans¬ 
mission  curve  C  becomes  tangent  to  the  h'f  curve  is  the  scaled  h7  at  (3) 
which  corresponds  to  the  same  height  (3)  on  the  oblique  ionogram  only 
shifted  by  sec  0.  The  curves  A  and  B  illustrate  how  further  points  on 
the  h  f  curve  may  be  obtained  by  shifting  the  transmission  curve. 

6.3.  The  Equivalence  Theorem 

Neglecting  the  earth' s  magnetic  field  and  applying  Snell' s  law,  the 
concept  of  equivalent  path  for  flat  earth  is  shown  in  figure  6.  3. 

The  equivalent  path  concept  states  that  the  time  required  for  a 
signal  to  travel  the  assumed  actual  path  TBR  is  the  same  as  a  wave 
traveling  at  the  speed  of  light  by  the  path  TAR  [Mitra  1952],  This  concept 
then  leads  to  an  important  relation  between  vertical  and  oblique  incident 
propagation;  i.e.  ,  the  height  of  the  equivalent  path  is  the  same  as  the 
vertical  height,  h/,  measured  at  the  equivalent  vertically  incident 
frequency. 

The  above  relationship,  secant  0  (corrected)  and  assumed  parabolic 
distributions  for  two  layers  (E  and  F2),  are  the  theoretical  basis  for  ray 
path  selections,  time  delays,  and  vertical  angles  of  arrival. 

6.  4.  Parabolic  Assumption 

The  parabolic -layer  method  used  is  based  on  the  equivalent  path 
concept  and  the  secant  law  for  curved  earth  and  ionosphere.  The  validity 
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Figure  ,.3.  Equivalence  Theorem  for  a  Plane  Ionosphere 

and  Plane  Earth 


of  this  method  is  determined  when  one  observes  how  well  the  shape  of  the 
ionosphere  is  approximated  by  two  parabolas.  It  is  believed  that  for 
average  conditions,  at  frequencies  of  interest  for  long-distance  HF 
propagation  and  considering  the  limiting  factor  of  accurate  world-wide 
data,  the  shapes  of  electron  density  versus  true  height  profiles  are 
adequately  represented. 

The  model  used  in  this  method  is  patterned  after  the  work  of 
Newbern  Smith  [  1939]  and  the  transmission  equations  developed  by  Rawer 
L1948,  1950],  [  Bibl  1950]  for  parabolic  ionospheric  layers.  The  effect 
of  the  earth*  s  magnetic  field  is  neglected  and  the  vertical  ionization 
distribution  is  assumed  to  be  parabolic.  The  earth  and  ionosphere  are 
assumed  concentric  and  the  vertical  ionization  distribution  taken  to  be 
constant  over  the  area  in  which  the  radio  ray  is  in  the  ionosphere.  The 
angle-of  elevation  at  the  transmitter  is  therefore  equal  to  the  angleof- 
arrival  at  the  receiver. 

The  range  equation  developed  by  Rawer  for  two  parabolic  layers  is 
(see  figure  6.  4) 

yFtan  aF  {  arc  coth  jaF 


Mrnar  0 


M-p 


-  1 


)] 


T  L  Rc 


/"arc  tanh  a. 
tan  ar  (  - ^ 


)]} 


(6.3) 


where 


.  sec  aE)|t  fE)F 
-  f 

D  =  grea*--circle  distance  -  km 

CCEjF  =  angle  of  incidence  at  the  midpoint  of  the  layer 

aE  r  -  arc  sin  |  —■  cos  A 
LRe>f  J 

Rn  =  earth*  s  radius  -  km 


A  =  angle  of  takeoff  -  degrees 
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f  =  operating  frequency  -  MHz 

i'E.F  -  critical  frequency  r  M'Hz 

^■E,F  “  Ho  +  hi-rx(£,F) 

^m{e  f)  =  true  height  of  maximum  ionization 

yE,F  -  semithickness  of  the  layer  in  kilometers. 

The  terms  of  the  range  equation  (6.  3)  correspond  to  the  angular 
distances  as: 

Term  1  is  0X  +  0n 

Term  2  is  0F 

Term  3  is  6a  (see  figure  6.4) 

The  above  range  equations  are  derived  for  one  hop  via  the  reflecting 
region.  In  the  case  of  multihop  propagation  the  method  applied  is  to 
convert  the  path  parameters  to  an  average  profile  along  the  path;  thus 
equal  hop  lengths  are  assumed  along  with  equal  angles  of  elevation.  Areas 
typical  of  actual  reflection  and  retardation  regions  are  inspected  to  gen- 
erate  the  typical  path  profile. 

Rawer' s  range  equation  and  the  transmission  curve  method  of 
Newbern  Smith  are  coupled  together  to  yield  a  fast  and  economical  method 
of  selecting  ray  paths  which  is  consistent  with  the  scaled  data. 

7.  Calculation  of  Sky-Wave  Paths 

The  calculation  of  the  upper  limit  of  frequency  and  signal  levels 
for  a  given  circuit  requires  knowledge  of  the  path  and  its  associated 
ionospheric  parameters  as  seen  in  the  preceding  sections.  The  total 
ground  range  (great-circle  distance)  of  the  circuit  is  calculated  by 
knowing  the  geographic  coordinates  of  the  transmitter  and  receiver.  The 
earth  is  considered  a  perfect  sphere  and  the  distance  is  calculated  by 

cos  (g)  =  sin  (xx  )  •  sm  (x2)  +  cos  (xx  )  (7.  0 

cos  (x2)  *  cos  (y2  -  yx) 
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where 


g  =  great-circle  distance  -  degrees 
xx  =  transmitter  latitude  -  degrees 
yx  =  transmitter  longitude  -  degrees 
x2  =  receivei  latitude  -  degrees 
y2  =  receiver  longitude  -  degrees. 

Once  this  great-circle  arc  is  defined  the  bearing  from  receiver 
to  transmitter,  degrees  east  of  north  (0  ^  b  ^  360  °),  is  calculated  by 

b  =  1 14.  5916  *  tan-1  •/[ sin  (u  -  90  +  X3)  •  sin  (u-g)  |  /  c  (7.  2) 

u  =  (180  -  xx  -  x2  +  g)/  2 
c  =  sin  (u)  •  sin  (u  -  90  +  ). 

With  the  great-circle  route  and  the  associated  bearings  defined, 
the  typical  reflection  and  penetration  areas  of  the  ray  in  the  ionosphere 
can  be  calculated  by 

Reflection  Lat  =*90  -  arccos  (x^)  (7.3) 

where 

x'i  =  cos  (P)  ‘  cos  (90  -  x2 )  +  sin  (p)  ♦  sin  (90  -  x3 )  •  cos  (b) 

p  =  great-circle  distance  in  degrees  from  transmitter 
terminal  to  reflection  latitude. 

The  geographic  longitude  of  the  reflection  area  is 

y  =  y2  -  y'  (7.4) 

and 

y'  =  arccos  j^cos  (p)  -  cos  (90  -  x2) 

•  cos  (90  -  x3) /  sin  (90  -  x2)  •  sin  (90  -  x3)  J 

Xjj  =  reflection  area  latitude. 
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The  zenith  angle  of  the  sun  at  any  typical  reflection  or  penetration 
area  on  the  great  circle  can  be  calculated  by 

cos  ('{')  =  sin  (2)  •  sin  (s)  +  cos  (z)  .  cos  (s)  •  (7.  5) 

cos  [15  X  GMT  -  180  -  y  J 

where 

^  =  sun' s  zenith  angle  -  degrees 

z  -  reflection  area  latitude  -  degrees  =  x3 

s  -  latitude  of  subsoiar  point  of  sun  for  middle 
of  month  in  question  -  degrees. 

The  above  parameters  are  needed  for  the  development  of  a  typical 
profile  of  electron  density  along  the  path. 

The  vertical  angles  of  arrival,  angles  incident  in  the  absorbing 
region,  and  angles  of  reflection  are  a  function  of  the  above  profile,  path 
length,  and  number  of  sky-wave  hops  involving  reflection  from  either  the 
E  or  F2  region. 

Seven  distinct  ray  paths  are  evaluated  for  each  hour  at  each  operating 
frequency.  The  ray  paths  must  be  geometrically  possible  with  a  takeoff 
angle  A  greater  than  any  arbitrary  values  set  by  the  user.  The  sky-wave 
paths  evaluated  are  (a)  two  ray  paths  reflected  by  the  regular  E-region 
(b)  three  ray  paths  reflected  by  the  F2  region,  and  (c)  two  mixed  modes 
of  the  "N"  type,  i.e.  ,  E  reflections  occurring  on  either  end  with  F2 
reflections  following. 

The  following  section  is  devoted  to  defining  the  angles  and  slant 
range  associated  with  the  above  ray  paths  using  the  transmission  curve 
method  and  the  parabolic  range  equations.  Figure  7.  i  illustrates  some 
typical  sky-v'ave  paths  which  are  evaluated  by  this  method. 
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7.  1.  Theoretical  Method  for  Computing  Virtual  Height  of  Reflection 
and  Angles  of  Arrival  for  Obliquely -Incident  Ray  Paths 
through  Parabolic  Layers 

Effective  communication  by  way  of  the  ionosphere  requires  a 
knowledge  of  the  range  of  frequencies  that  can  be  propagated  between  two 
points  on  the  earth' s  surface. 

A  knowledge  of  the  elevation  angle  at  each  frequency  is  also  impor¬ 
tant  when  directive  antennas  are  employed  to  permit  reduced  transmitter 
power.  A  transmission  curve  is  "a  simple,  rapid  graphical  method  of 
obtaining  from  vertical-incidence  data,  he  maximum  usable  frequency 
over  a  given  path,  and  the  effective  heights  of  reflection  (hence  elevation 
angles)  of  waves  incident  obliquely  upon  the  ionosphere.  " 

The  quotation  above  is  from  a  1939  Proceedings  of  the  I. R.E.  paper 
in  which  Newbern  Smith  presented  the  geometrical  basis  and  the  mathe¬ 
matics  of  transmission  curves  for  various  combinations  of  flat  and  curved 
earth  and  ionosphere.  The  paper  includes  references  to  the  earlier  work 
of  Breit  and  Tuve,  Martyn,  Eckersley,  Berkner,  Millington,  and  others. 
Transmission  curves  ace  also  discussed  in  'Ionospheric  Radio  Propa- 
tion"  NBS  Monograph  80.  (Davies  1965,  Chapter  4]. 

This  section  shows  how  one  of  Smith' s  many  equations  may  be  rear¬ 
ranged  to  provide  a  rapidly  convergent  method  which  imitates  the  appli¬ 
cation  of  a  transmission  curve  to  a  vertical  incidence  ionogram.  At  present 
the  ionogram  is  assumed  to  be  adequately  represented  by  a  parabolic 
distribution  of  electron  density,  but  other  models  for  empirical  distributions 
could  be  incorporated  in  the  routine.  Several  model  distributions  are 
considered  in  section  3.  3.4  of  'Ionospheric  Radio  Propagation," 

[  Davies  1965]. 

Except  for  notation.  Smith' s  equation  giving  an  approximate*  equiv¬ 
alence  relation  for  curved  earth  and  ionosphere  is  as  follows: 

*  The  approximation  is  that  both  h'  -h^,  and  h-h0  be  less  than  400  km  and 
hence  much  less  than  R0  +  h0  . 
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tan2  0 


(7.6) 


f’  -  f  cos  0  \T~~  2  (h1  -h) 

'  R0  +  ho 


where  {see  figure  7.2): 

i'  =  the  equivalent  vertical-incidence  frequency,  i.  e.  ,  th§ 
frequency  on  the  vertical  incidence  ionogram  which 
correponds  to  f 

f  =  the  probing  frequency  at  oblique  incidence 

0  =  the  half -vertex  angle  of  the  equivalent  triangular  path 
composed  of  the  straight  lines  joining  the  transmitter 
and  receiver  locations  with  a  point  at  the  height  h' 
above  the  midpoint 

h'  =  the  virtual  height  on  the  vertical  incidence  ionogram 

at  the  frequency  f  ,  measured  above  the  earth' s  surface 

h  =  the  true  height  of  reflection  of  the  curved  ray  path  in  the 
ionosphere 

hQ  =  the  height  of  the  bottom  of  the  ionospheric  layer 

R0  =  the  earth' s  radius 

D  =  the  surface  distance  between  the  end  points  of  the  hop 

Equation  0.6)  is  basically  the  well-known  "secant  law"  which  holds 
for  a  flat  earth  and  ionosphere.  The  square  root  represents  k,  a 
"correction  factor"  resulting  from  the  curved  geometry.  This  equation 
is  sometimes  written  as 


f  =  k  f  sec  0  ,  and  (7  7) 

comparison  with  equation  (1)  shows  that  this  correction  factor  is 


1/k 


^/7  Z{h'  -h) 
'  R0  +  K 


tan2  0  . 


(7.8) 


The  value  of  k  ranges  from  unity  for  very  short  distances  to  about  1.  2 
for  the  maximum  one -hop  distance  via  the  F  layer  of  about  5000  km. 
Since  k  is  a  function  of  0,  h„  ,  and  h'  -h,  it  depends  on  the  distance  of 
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transmission  and  the  height  and  the  electron  density  distribution  of  the 
ionosphere.  In  practice  an  approximate  correction  factor  which  depends 
only  on  the  distance  (figure  6.  1)  in  used  to  simplify  the  construction  of 
the  transmission  curves.  A  logarithmic  transmission  curve  is  illustrated 
in  figure  7.3  taken  from  Davies  [1965]  which  also  describes  the  derivation 
and  application  of  these  transmission  curves.  Kobayashi  [1961]  has 
developed  transmission  curves  using  the  explicit  dependence  of  k  on  h'  -h 
and  showing  how  an  ionogram  may  be  corrected  to  allow  for  the  effect  of 
ionospheric  curvature.  The  method  described  here  also  allows  for  the 
explicit  dependence  of  k  on  h-h'  . 

The  iterative  formula  for  the  computer  routine  is  obtained  from 
equation  (7.  6)  as  follows. 

Given  are  the  oblique  incidence  frequency,  the  distance  betv/een  the 
end  points  of  the  hop  and  the  critical  frequency,  maximum  layer  height 
and  semithickness  at  the  midpoint  (M  in  figure  7.  2).  These  parameters 
are  sufficient  to  specify  the  midpoint  vertical-incidence  ionogram  mathe¬ 
matically  through  parabolic  layer  theory.  For  the  parabolic  layerf  the 
virtual  height  is  given  by 

h'  =  h0  +  y  x  arctanh  (x)  [Appleton  and  Beynon  1947],  (7.9) 

where 


y  =  hTn  -  h0  =  the  semithickness 


h^,  =  height  of  maximum  electron  density  of  the  layer 


x 


the  vertical!  incidence  sounding  frequency  normalized 
to  the  critical  frequency 


fc  =  vertical  incidence  critical  frequency  (corresponds  to 
maximum  electron  density). 
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for  Curved  Ionosphere 


(7.11) 


The  true  height  of  reflection  in  a  parabolic  layer  is 

h  =  h0  +  y  (1  -  /  1  _  x2  )  >  (7.  10) 

and  the  difference  between  the  virtual  and  true  height  is 

h*  -h  =  y  (  n  -  x2  +  x  arctanh  x  -1)  =  Y  5.  (7. 

Both  the  above  expressions  refer  to  vertical  incidence  at  the 

f ' 

midpoint  and  involve  x  =  ,  the  quantity  sought. 

4 

f'  /f'\  /f  N 

Now  write  —  =  \^£~J  normalizes  both  f'  and  f  to  the  mid¬ 

point  critical  frequency.  Rewriting  equation  (7.  6)  yields 


(DC 


4a  * 

— •  )  =  COS  0 


1  -  -  tan2  0 

R0  +  ^ 


(7.12) 


Multiplying  both  sides  by  —  and  then  squaring  both  sidjs  yields 

^ADDDD1-^  tan2  0  ^  cos2  (7.13, 

For  convenience  in  computing  we  express  cos  0  in  terms  ox  tan 

0:  cos2  0  -  - - — 5— —  .  From  the  geometry  of  figure  7.2,  it  is  found 

1  +  tarn  0 


tan  0  = 


Rn  sin  0 


sin  0 


R0  (1  -  cos  0)  +  h'  1  -  cos  0  +  h1 


(7.14) 


where  9  -  — —  ,  the  central  angle  between  the  midpoint  and  either  end 
2R0 

of  the  hop. 

In  drawing  figure  7.  2  with  the  0  vertex  closed,  it  is  assumed  that 
the  height  of  the  equivalent  triangular  path  xs  equal  to  the  virtual  height 
h'  corresponding  to  f'  on  the  ionogram.  This  is  true  for  the  flat  iono¬ 
sphere  [Waterman  1952].  However,  in  this  practical  work,  the  error 
due  to  this  difference  is  a  second  order  effect. 


50 


Equation  (7.  14)  for  tan  0  involves  the  distance  of  transmission  and 
the  virtual  height  which  is  given  by  equation  7.  9.  The  complete  equation 


to  be  iterated  is 
2 


2 

x  r  = 


CrJi.  CO  0- IS 


1  +  tan2  0 


1  + 


c- 


sin  9 


1  -  cos  9  + 


ho  +  y  Xo  arctanh  Xo 


-)! 


(7.15) 


where 

6  =  X,.  arctanh  x.-l+Zl-Xe2  {from  7.  1 1)  ;7.  16) 


and 


xa  =  an  estimate  of  x 

xr  =  the  result  of  computing  the  r.h.  s.  of  (7.  15)  using  x<, 
(Note  that  x,  depends  mainly  on  0  and  h^). 


To  start  the  iterative  procedure  we  use  the  limits  x,  =  0  and 
x8  =  xt  ,  where  xt  represents  the  point  on  o\ir  mathematical  h'  -  f  iono- 
gram  tangent  to  a  transmission  curve  positioned  for  finding  the  maximum 
usable  frequency.  This  upper  limit  also  assures  that  only  values  of  h' 
corresponding  to  low  angle  rays  are  found  by  this  method. 

With  the  lower  limit 


Xo  arctanh  Xe  =0,  6=0  (hence  k  =  1)  and  equation  (7.  15}  reduces  to 


sin  0  g  (7.17) 

1  -  cos  0  +  ~  J 
1\0 


which  depends  only  on  the  transmitter  distance  and  the  height  of  the  bottom 
of  the  parabolic  layer.  The  resulting  xr  is  already  closer  to  but  less  than 
the  true  value  which  results  from  further  iterations. 
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From  some  other  work  [Remmler  1965]  it  is  found  that  the  upper 
limit,  depends  only  on  y/  h*,  as  follows: 


(7.18) 


Thus  xt  decreases  as  increases.  For  the  F  layer,  y/ h0  varies 

fro 

between  about  0.  2  and  1. 0;  hence  xt  varies  between  about  0.  92  and  0.  72, 
and  the  corresponding  6  varies  between  0.854  and  0.347,  respectively. 
Therefore,  for  any  reasonable  value  of  y,  the  subtractive  term  in  the 
numerator  varies  less  than  one -tenth  as  fast  as  the  additive  term  in  the 
denominator  of  equation  (7.  15).  Alternatively,  ~  varies  between  1  and 
0.  7  as  k  varies  between  its  limits  of  1  and  1.2.  Finally,  since  tan  0 
varies  from  zero  at  vertical  incidence  to  about  3  at  the  one -hop  F  limit 
the  entire  denominator  can  vary  from  1  to  about  11.  Hence,  xr  is  deter¬ 
mined  principally  by  the  transmission  distance  and  layer  height  as  indicated 
above . 

If  the  upper  limit  xt  is  used  as  an  estimate  of  x  in  (7.15)  the 
resulting  value  of  xr  is  closer  to  but  greater  than  the  true  value.  A 
second  iteration  using  either  of  the  above  values  of  x^  gives  a  better 
approximation  which  always  lies  on  the  same  side  of  the  true  value  as 
the  input  estimate.  Faster  convergence  occurs  when  the  next  estimate 
Xo  is  the  average  of  the  first  two  xr'  s.  Further  estimates  are  obtained 
by  averaging  each  new  result  with  the  closer  of  the  two  values  just  pre¬ 
viously  averaged.  This  procedure  gu'.rantees  that  the  two  values  averaged 
straddle  the  unknown  true  value.  Virtual  heights  are  computed  for  each 
x,  and  when  the  difference  in  successive  virtual  heights  is  less  than  0.  5 
km  the  process  is  stopped.  The  average  virtual  height  is  then  used  to 
find  the  angle  of  elevation  of  the  F-region  ray. 
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Three  or  four  iterations  are  usually  sufficient  to  reduce  the 
virtual  height  difference  to  0.  5  km.  If  xr  exceeds  xt  when  xt  is  put  into 
the  r.h.  s.  of  (7.  15)  the  skip  distance  exceeds  the  transmitter -receiver 
distance,  that  is  f  exceeds  the  MUF  for  the  path. 

7.  2.  Method  for  Incorporating  the  Bending  in  the  E  Layer 
in  Predicting  Ray  Faths  by  V/ay  of  the  F2  Layer 

Rawer  [i960]  and  Woyk  [  1959]  among  others  have  discussed  the 
refraction  (or  bending)  of  radio  waves  passing  completely  through  a 
spherical  ionized  layer  using  a  simple  ray  treatment.  Figure  7.4 
illustrates  the  ray  geometry  and  the  nomenclature  employed. 

The  following  formulas  and  data  are  supposedly  given  (at  this  point): 

1.  Rawer' s  formula  for  the  bending  3  experienced  by  a  ray 
passing  entirely  through  a  parabolic  layer: 

p  =  2ctana(tiSlpi),  (7.19) 

where 

Ye 

e  =  F —  u — pr  :  yE  “  semithickness  of  E  layer 
+  hmE 

h^E  =  height  of  maximum  density 
of  E  layer 

R0  =  radius  of  earth 

a  =  angle  of  incidence  of  the  unrefracted  ray 
extension  at  hj^E  measured  from  vertical 

fE 

u  =  —  sec  a 

fE  =  critical  frequency  of  E  layer 

f  =  signal  frequency. 

Rawer  and  Woyk  have  shown  that  3  can  also  be  measured  at  the 
center  of  the  earth,  i.  e.  ,  aetween  perpendiculars  to  extensions  of  the 
rays  above  and  below  the  layer  penetrated  (figure  7.4). 
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A  =  ELEVATION  >'  V-iLE  AT  T*  (NO  BEND'NG) 

A/3  -  ELEVATION  »N6LE  AT-Tx  (WITH  BENDING) 
a  =  ANGLE  OF  INCIDENCE  AT  hm  E  (NO  BENDING) 
ap  -  ANGLE  OF  INCIDENCE  AT  hm  E  (WITH  6ENOING) 

4>  :  ANGLE  OF  INCIDENCE  AT  h’F2  (NO  B1NOING  ! 

4>P  =  ANGLE  OF  INCIDENCE  AT  h’  F2  (WITH  BENDING) 

/?  ■  BENDING  IN  E-REGION 

R)on--  RADIUS  OF  TANGENCY  OF  BENT  RAY  EXTENSIONS 
R0  =  RADIUS  OF  EARTH;  AIL  HEIGHTS  ( h  ,  h‘  )  ABOVE  EARTH’S  SURFACE 
Tx  =  TRANSMITTER 
D/2  =  DISTANCE  BFTWEE'I  Tx  AND  M 
$  ■■  O/ZRo  ,  er-d-P 

Figure  7.4.  Geometry  of  E-Region  Bending  for  Rays 
Propagating  by  the  F2  Layer 
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2.  The  iterative  routine  finding  the  virtual  height  of  reflection 

in  the  F  layer  as  a  function  of  the  F-layer  parameters  and  one 
hop  path  length.  This  routine  is  described  in  detail  in  section 
7.  1. 

3.  Snell' s  law  for  spherical  layers  [Woyk  1959;  Rawer  I960  ]; 

n  p  sin  a  =  const.  ,  (7.20) 

where  n  -  index  of  refraction  and  p  =  radius  of  curvature 
where  a  is  measured  w.r.t.  the  vertical  at  that  radius. 


Whenever,  n  =  1,  e.g.  ,  between  layers,  we  have 

p  sin  cc  =  const. 

This  applies  to  straight  line  extensions  into  ionized  regions  also.  The 
constant  may  be  determined  from  any  corresponding  pair  of  values  p  and 
a,  or  R0  +  h'  and  0,  etc.  Thereafter  a  can  be  found  for  any  given  p  and 
vice  versa. 

4.  The  following  parameters  for  both  E  and  F2  layers: 

a.  critical  frequency,  fE ,  fF 

b.  height  of  bottom,  hoE,  h^,  F 

c.  semithickness,  yE,  yf 

d.  height  of  max  density,  h^,  =  ho  +  v. 


5.  The  signal  frequency  (f)  and  the  hop  length  (D),  expressed  in 
radians  and  measured  at  the  center  of  the  earth,  i.  e.  ,  : 

§  =  R„e  (7.21) 


or 


\ 

a 


D 


the  central  angle  from  either 
endpoint  to  the  midpoint. 
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The  iteration  proceeds  as  follows: 

A.  Use  the  virtual  height  routine  (2  above)  with  0,  f  and  the  F  layer 
parameters  to  determine  the  elevation  angle  (A)  for  propagation  by  the 
F  layer  alone.  This  angle  can  be  found,  for  example,  through  application 
of  (7.20)  above  as  follows: 

R0  sin  (~  -  A)  =  (R0  +  h'F)  sin  0  (7.22) 

or 

R0  cos  A  =  (R0  +  h'F)  sin  0 
A  Ro  +  h'  f  ^ 

cos  A  =  - - - -  sin  0. 

Ro 

R.  The  angle  of  incidence  (a)  at  h^  of  the  i.  layer  is  obtained  from 
(7.  20)  in  the  same  manner: 

(R0  +  hmE)  sin  a  =  Rc  sin  (tt/  2  -  A) 


fR0  sin  (n/  2  -  A)^ 
a  =  arcsm  (  — - - — — — -  ) 

^  R0  +  h^E  J 


(7.23) 


or 


C.  Next  a  from  step  3  and  the  E  layer  parameters  fe ,  h^E,  and 
yE  are  substituted  into  Rawer' s  formula  (7.  19)  to  get  an  initial  value  of 
the  bending,  3. 

D.  The  ray  which  would  reach  the  receiver  by  way  of  the  F  layer 
alone,  is  bent  in  passing  through  the  E  layer;  it  overshoots  the  receiver 
location.  Now  a  ray  leaving  the  transmitter  at  a  somewhat  greater 
elevation  angle  can  be  found  for  which  the  .angular  distance  added  by  the 
bending  is  just  compensated  by  the  reduction  in  angular  distance  below 
and  above  the  E  layer  caused  by  increasing  the  elevation  angle  (figure  7.  4). 
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Because  the  bending  is  itself  a  function  of  the  elevation  angle  as  is  the 
virtual  height  of  reflection  in  the  F  layer,  a  successive  approximation 
procedure  is  used.  Steps  A,  B,  and  C  initiated  the  process  and  similar 
steps  are  used  in  the  iteration. 

E.  Subtract  the  total  initial  bending  from  9  to  find  the  reduced  half¬ 
hop  distance: 

6r  =0-9.  (7.24) 

Repeat  steps  A,  B,  C  with  9P  substituted  for  0  to  find  a  new  smaller 
bending  corresponding  to  the  decreased  angle  of  incidence  for  the  reduced 
distance.  The  subscript  Pin  figure  7.4  indicates  angles  with  some 
bending. 

F.  Since  this  bending  is  less  than  anticipated,  the  ray  falls  short 
of  the  receiver  and  the  final  bending  lies  between  the  two  values  computed 
thus  far.  Therefore,  average  these  two  values  of  bending,  compute  the 
corresponding  reduced  distance  as  in  step  E  and  iterate  through  steps  A, 

B,  C,  and  F  until  the  F  layer  virtual  height  of  reflection  falls  within  the 
prescribed  limits,  (or  possibly  until  the  sum  of  the  reduced  distance  and 
the  bending  is  within  a  prescribed  small  angular  increment  of  the 
desired  distance). 

8.  Transmission  Losses 
8.  1.  Absorption  in  Lower  Layers  (Lj) 

The  main  features  of  the  structure  and  the  diurnal,  seasonal,  and 
geographical  variations  of  the  D  and  E  regions  are  adequately  explained 
in  terms  of  the  ionization  produced  by  absorption  of  ultraviolet  solar 
radiation  and  the  recombination  of  the  ions,  and  are  assumed  to  be  related 
to  the  zenith  angle  of  the  sun.  The  absorption  equation  of  section  8.  2  is  a 
semiempirical  relationship  derived  from  observations  of  operating  high 
frequency  circuits. 


The  mean  effective  height  of  the  absorbing  region  is  taken  to  be 
100  km  above  the  earth  and  the  total  absorption  of  a  high  frequency  wave 
is  directly  related  to  the  length  of  the  trajectory  within  the  absorbing 
region  [Laitinen  and  Haydon  1962]. 

In  this  development  of  a  semiempirical  expression  for  absorption 
in  the  lower  layers,  it  was  assumed  that  there  is  negligible  nighttime 
absorption  in  the  range  of  frequencies  3  to  30  MHz  so  that  the  change 
in  signal  level  throughout  the  day  was  the  basis  for  the  analysis.  Likely 
propagation  paths  in  the  daytime  were  estimated  from  median  ionospheric 
conditions  without  any  adjustment  for  the  likelihood  of  the  various  sky-wave 
paths.  Comparison  between  circuit  performance  predictions  based  on 
monthly  median  ionospheric  conditions  and  the  observed  circuit  perform¬ 
ance  indicated  that  the  median  signal  level  was  8.  9  dB  below  the  quasi¬ 
maximum. 

The  prediction  method  presented  in  this  report  considers  the  prob¬ 
ability  of  a  sky-wave  path  and  combines  this  probability  with  the  expected 
day-to-day  signal  level  variations  during  the  days  a  sky-wave  path  is 
expected.  Since  the  signal  level  distributions  required  in  this  report  apply 
only  to  those  days  a  sky-wave  path  is  expected,  the  adjustment  between 
the  median  calculated  for  these  days,  and  the  expected  signal  level  should 
be  somewhat  less  than  previously  required.  Pending  the  completion  of 
more  detailed  comparisons,  the  8.  9  dB  is  retained  for  the  temperate 
regions.  The  values  for  other  latitudes  are  shown  in  Appendix  C. 

8.2.  Theoretical/ Empirical  Determination  of  Lower  Region  Absorption 

The  absorption  that  takes  place  in  the  lower  E  and  upper  D  region 
(80-i20  Jem)  of  the  ionosphere  is  assumed  to  account  for  much  of  the  loss 
encountered  by  a  ray  traveling  between  transmitter  and  receiver  via  the 
ionosphere.  Nondeviative  absorption  (u  -  1)  is  the  only  absorption 
explicitly  dealt  with  'n  these  calculations  of  monthly  median  transmission 


58 


losses.  Deviative  absorption  is  in  part  considered  in  the  distributions 
(Yp)  about  these  monthly  median  values.  The  expression  used  for  the 
absorption  is  semiempirical  and  a  brief  explanation  of  its  origin  is  given. 
A  complete  explanation  of  data  and  analyses  is  contained  in  Technical 
Report  No.  9  [Laitinen  and  Haydon  1962].  Appleton  derived  an  expression 
for  the  total  nondeviative  lower  region  absorption  for  frequencies  much 
greater  than  the  collision  frequency.  The  equation  expressed  in  terms  of 
the  ionospheric  reflection  coefficient  for  vertical  incidence  is 


-  loge  P  =  4.  13H 


4n2e2Nr 


cos 


1.5 


mc(w  +  )2 


(8.1) 


where  p  =  the  reflection  coefficient 
e  -  the  electronic  charge 
m  =  the  mass  of  the  electron 
c  =•  the  velocity  of  light 

N0  =  the  electron  density  at  the  maximum  ionization  of  the  layer 
when  the  sun' s  zenith  angle  is  zero 

vQ  =  the  collisional  frequency  at  the  height  of  maximum  ionization 

H  =  the  scale  height 

'll  =  the  zenith  angle  of  the  sun 

=  the  angular  frequency  of  the  propagated  radio  wave 

i»L  the  angular  gyro -magnetic  frequency  due  to  the  longitudinal 
component  of  the  earth's  magnetic  field. 


This  equation,  based  on  a  Chapman  distribution  of  ionization,  shows  the 
vertical  incidence  absorption  for  the  ordinary  component  of  the  electro¬ 
magnetic  wave.  The  absorption  varies  directly  with  cos  '*5  and 
inversely  with  the  square  of  the  sum  of  the  wave  frequency  and  the  longi¬ 
tudinal  component  of  the  gyro -frequency. 

The  above  equation  indicates  that  absorption  would  cease  at  a  solar 
zenith  angle  of  90°.  The  data  analyzed  indicated  that  it  in  fact  levels  off 
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at  sunset  at  the  ionospheric  height  of  .110  kin  or  a  zenith  angle  of  102 
The  portion  of  the  equation  is  then  adjusted,  '.o  be 

(cos  K  <ii)m  =  (cos  ()„  881  t)m  (8.  2) 

K  =  constant  of  90°/  102.  2 0  =  0.  881 

m  =  constant  obtained  from  data  =1.3 
=  zenith  angle  of  the  sun. 

Analyses  of  the  data  over  a  period  of  several  years  indicated  that 
an  adjustment  needs  to  be  made  to  allow  for  the  change  in  solar  activity. 
Twelve-month  running  average  sunspot  numbers  were  used  to  analyze 
the  data  in  an  attempt  to  remove  the  seasonal  trends  of  absorption. 

The  analysis  produced  a  linear  relationship  with  sunspot  number: 

As  =  A,  (1  +  b  S)  ,  (8.3) 


where 

S  =  Twelve -month  running  average  sunspot  number 
b  =  constant  of  0.  0037  obtained  from  the  experimental  data 
Aq  =  absorption  at  sunspot  number  zero 
As  =  absorption  at  sunspot  S- 


The  total  equation  for  the  anticipated  losses  (A)  due  to  the  absorption 
in  the  lower  regions  on  oblique  paths  is  written  as 


A  = 


615.  5  n  sec  g>[  1  4-  0.  0037S]  [cos  0.  881  fj1 *3 ° 
[f  +  fn]1’38 


,  (dB) 


(8.4) 


where 

sec  0  =  factor  relating  back  to  vertically  incident  waves 

n  =  number  of  hops  in  the  ray  path  from  transmitter  to 
receiver 

fH  =  gyro-frequency  at  100  km  above  earth’ s  surface  -  MHz. 
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The  above  equation  has  the  same  form  as  the  Appleton  equation,  but 
contains  the  aforementioned  correction  in  addition  to  the  factor  615.  5 
determined  from  oblique  observations  and  the  fH  used  in  conjunction  with 
the  modified  exponent  of  1.98.  Errors  in  the  above  equation  can  develop 
(a)  because  the  frequency  is  not  sufficiently  greater  than  the  collision 
frequency  or  (b)  at  those  times  when  the  (f  +  fH)1,98  expression  does  not 
sufficiently  approximate  (f  +  ft)2  where  it  is  the  longitudinal  component  of 
the  magnetic  field  for  the  path  in  question. 

It  has  also  been  stated  that  the  variation  of  absorption  between 
winter  and  summer  is  less  than  that  expected  by  theory  based  on  solar 
zenith  angle.  It  was  found  by  the  data  anal) ^ed  that  the  absorption  follows 
an  inverse  relationship  with  layer  height.  During  the  summer,  layer 
heights  are  typically  higher  than  in  the  winter;  thus  the  ray  angles  at  the 
absorbing  regions  are  higher  for  a  given  ray  path.  This  in  turn  tends  to 
decrease  absorption  during  the  summer  months  relative  to  the  values  during 
the  winter  months;  thus  explaining  the  difference  between  seasons  not 
accountable  explicitly  by  the  solar  zenith  angle  dependence.  Equation  8.4 
is  a  reasonable  fit  to  the  data  of  Technical  Report  No.  9  [.Laitinen  and 
Haydon  1962]  for  values  of  f  +  fH  in  excess  of  5  MHz;  however,  the  equation 
fails  to  fit  the  data  at  lower  values  of  f  +  fH  .  A  revised  equation  which 
yields  similar  results  as  equation  8.4  for  values  of  f  +  fH  greater  than 
5  MHz  and  also  represents  the  data  to  a  lower  limit  of  3  MHz  is 


.  t77.  2  n  (sec  <P)  I 

dE  ~  ff  +  f„)1.88  +  10.2 


(8.  5) 


where 


I  =  [l  +  0.00373]  [cos  .881  f]1-30. 
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The  analysis  of  observations  of  nighttime  field  intensities  has 
indicated  that  the  lower  region  absorption  does  not  in  fact  cease,  but  rather 
tends  to  level  off  as  the  absorption  index  (defined  p.  8)  I  approaches  a 
value  of  one -tenth.  The  above  absorption  equation  is  thus  restricted  to 
values  of  I  greater  than  0.1.  This  assumption  checks  well  with  the  work 
of  Wakai  [1961]  and  unpublished  comparisons  with  observations  at  The  In¬ 
stitute  for  Telecommunication  Sciences  and  Aeronomy  (formerly  CRPL). 

For  a  complete  explanation  of  the  above  analysis  and  data,  refer  to 
Technical  Report  No.  9  of  the  U.  S.  Army  Radio  Propagation  Agency 
[Laitinen  and  Haydon  1962]. 

8.3.  Losses  Due  to  ground  Reflections  (L^) 

The  ground  reflection  losses  calculated  in  this  method  are  for 
randomly  polarized  sky  waves  which  are  assumed  to  have  equal  amounts 
of  energy  in  the  horizontally  and  vertically  polarized  fields.  The  losses 
are  represented  by  the  following  equation: 

k  =  10  logjo  Kit*!  (dB).  (8.6) 


v’here 

Kv  =  magnitude  of  the  vertical  reflection  coefficient  at  the  angle  (A) 

Kh  =  magnitude  of  the  horizontal  reflection  coefficient  at  the 
angle  (A) 

Kv  and  KH  equations  can  be  found  in  section  8.  5. 

It  should  also  be  mentioned  here  that  "cutback  factors"  [Wait  and  T.onda  1958] 
affecting  the  radiation  at  low  angles  are  not  at  the  present  time  included 
in  the  ground  reflection  losses,  and  it  is  assumed  that  the  majority  of  the 
energy  is  reflected  at  an  angle  of  reflection  equal  to  the  angle  of  incidence 
of  the  wave;  plans  are  under  way  for  improving  this  aspect  of  the 
predictions. 
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8.4.  Free  Space  Loss  Between  Isotropic  Radiators  (Lbf) 

Consider  isotropic  radiators  in  free  space  which  are  radiating  a 
power  w  (watts)  which  produce  a  field  intensity  of  w/  4trd2  when  d  is  much 
greater  than  X.  The  absorbing  area  (a)  of  an  isotropic  radiator  in  free 
space  is 

a  =  Xs  /  4 tt  ,  (8.  7) 


with 


X  =  free  space  wavelength. 


The  resulting  power  available  from  an  isotropic  radiator  at  a 
distance  d  from  the  isotropic  transmitting  antenna  is 


w,  =  w 


(4) 


(8.8) 


and  the  resulting  basic  transmission  [Norton  1959]  between  isotropic 
radiators  can  be  expressed  by 


=  10  logiQ  Qrr/  w.  ^ 


or 


(8.9) 


=  32.  45+  20  log1Qd  +  20  log^f 
f  =  frequency  -  MHz 

d  =  distance  in  kilometers  along  the  ray  path  between  the  two 
antennas. 


8.  5.  Theoretical  Antenna  Power  Gain  Equations  (GT  ,  GR) 
Equations  are  given  for  the  power  gain  relative  to  a  free  space 
isotropic  source  for  several  antennas  over  finite  earth  in  Appendix  A. 


*  Some  symbols  appearing  in  this  section  and  Appendix  A  may  duplicate 
symbols  used  for  different  variables  appearing  in  other  sections  of  the 
report.  It  was  necessary  to  do  this  when  attempting  to  remain  consistent 
with  the  symbols  used  in  antenna  theory. 
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This  section  contains  a  brief  presentation  of  ground  reflection  coefficients, 
power  gain  definition,  the  radiation  vector  method  and  ground  reflection 


factor. 


8.  5.  1.  Reflection  Coefficients 


If  an  electromagnetic  wave  in  air  is  incident  on  the  ground,  then 
the  reflection  coefficient  of  the  wave  is  defined  as  the  ratio  of  the  field 
strength  of  the  reflected  wave  to  the  field  strength  of  the  incident  wave. 
Since  the  incident  and  reflected  waves  have  phase  as  well  as  amplitude, 
the  reflection  coefficient  is,  in  general,  a  complex  number.  For  waves 
polarized  parallel  to  the  plane  of  incidence  of  the  wave  (vertical  polar¬ 
ization)  the  reflection  coefficient  is  [Schelkunoff  and  Friis  1952] 


Kv  =  - 


_  (er  -  ix)  sin  A  -  / (er  -  ix)  -  cos^  A 


(8.10) 


(Sr  -  ix)  sin  A  +  /(c^  -  ix)  -  cos2  A 


and  for  waves  polarized  perpendicular  to  the  plane  of  incidence  (horizontal 
polarization) 


K  =  sin  A  -  /(er  _  ix)  _  cos2  A 
sin  A  +  /(er  -  ix)  -  cos2  A 


(8.11  ) 


where  /  denotes  the  principal  branch  of  the  complex  square  root  function, 


ef  =  relative  dielectric  constant  of  earth, 
x  =— —  =  18  x  103  a/f 


=  dielectric  constant  of  free  space  (farad/ meter) 
V 

a  =  conductivity  of  earth  (mhos/  meter) 

U)  =  angular  frequency 

f  =  frequency  in  megahertz 
A  =  angle  of  elevation  in  degrees,  and 
i  =/-l  . 


If,  using  the  appropriate  subscript  V  or  H  on  K,  one  writes  K  in 
the  form  K  =  |KJ  elY,  then  JK|  is  the  amplitude  of  the  reflection  coefficient 
and  y  the  phase .  Since  the  function  elY  is  periodic  with  2TTi,  the  expression 
given  here  does  not  uniquely  define  the  phase.  It  will,  however,  always 
be  clear  from  the  context  which  value  of  Y  is  being  used.  Letting 
*  =  Y  ±  tt  , 

K  =  -  lKle1*  .  (8.12) 


It  is  the  purpose  of  this  section  to  obtain  K  and  >|ifor  both  vertically  and 
horizontally  polarized  waves.  Letting  A  =  (ep  -  ix)  -  cos2  A,  one  finds 
that 

ia 

A  =  pe  ,  (8.13  ) 

where 

JL 

P  =  C(ep  -  cos3  A)2  +  x3]2  ,  (8.14) 


and 


a  =  -Tan"1 


x _ 

(ep  -  cos2  A) 


(8.15) 


Equation  (8.  11)  can  now  be  written  as 

x. 

(ep  -  ix)  sinA  -  p2  e 

Kv  - - x 

(s  p  -  ix)  sinA+  pa  e 


i  a/  2 


i  a/  2 


(8.16) 


The  reduction  of  (8.16)  to  the  form  (8. 12)  gives 


J^p2  +  +  x2)2  sin4  A  -  2p  (ep2+  x2)  sin2  Acos  ( 

'a  +  ?.  Gin"1  x  ^l2 

v.  is  J 

(6V2  +  x2)2  (8.  !7) 

r  -  - 

[p  +  (sr2  +  x2)  sin2  A  +  2p2  (ep2  +  x3)2  sin  Acos  ( 

vj  +  sin  ,  » 

(e,3  +  x3  )2 
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The  analysis  of  t|fr  yields  the  following  three  cases: 


1.  If  (ers  +  x3)A  -  p  =  0,  then 

a.  If  sin  A  sin  (Tan-1  ~  +  r-)  -  0,  then  ^  =  0, 

r  ^ 

b.  If  sin  A  sin  (Tan-1  —  +  —  )  >  0,  then  ~  . 

c.  If  sin  A  sin  (Tan-1  —  +~)  <0,  then  %  =  -  ~  . 

Ep  u  L 


(8.18a) 

(8.18b) 


(8.18c  ) 


2.  If  p  -  (€r2  +  x2)  sin2  A  >  0,  then 


,  rZP2  /s~2  +  x*5  sin  A  sin  (tan-1  “  +  %  )  “] 
i|»v  =  Tan  1  _ r  _  er  2  (8.18d  ) 

-  p  -  (era  +  x2)  sin2  A  J 


3.  If  p  -  (ef2  +  x2)  r-in2  A  <  0,  then 


>{rv  -  Tan 


■i*  r* 

j-2?2  /Aet2  +  xs  sin  A  sin  (Tan-1  —  ) 

_1  j _ £r  L 

P  -  (f/  +  x2)  sir2  A 


+  n 


(8.18e  ) 


A  similar  analysis  leads  to 


1 

[  p*  -h  sin4  A  -  2p  sin2  A  cos  a  ]  2 

[p+  sin2  A  +  2p2  sin  A  cos  Hi] 

2 


(8.19  ) 


♦h  =  Tan-1  A  SL 2  ] 

L  P  -  sin2  A  ! 

1 —  _/ 


(8.20  ) 
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8.  5.  ?..  Power  Gain  Definition 


Following  Ramo  and  Whinnery  [i960],  the  gain,  g,  relative  to  an 
isotropic  source  is  defined  as  the  ratio  of  power  required  from  the 
isotropic  source  to  produce  the  given  intensity  in  the  desired  direction 
to  that  required  from  the  actual  antenna: 

4TTr2  Pr  4ttr2  P, 

S  “  w  "  iRP  *  (8.21  ) 


where 

Pr  =  the  time  average  Poynting  vector, 

r  =  the  distance  from  the  origin  of  the  coordinate  system 
to  point  in  space, 

w  =  the  total  radiated  power  in  free  space, 

R  x  the  radiation  resistance  of  the  antennr.. 

J.  =  the  input  current  to  the  antenna. 


8.  5.  3.  Radiation  Vector  Method  and  Ground  Reflection  Factor 

Using  the  radiation  vector  approach  introduced  by  Schelkunoff 
[Ramo  and  Whinnery  I960], 

pr  =  51F75  {  |Ne!2  +  |N0|S  ]  ,  (8.22) 

where 


X  =  the  wave  length  in  meters, 

T|5:120  tt  =  the  intrinsic  impedance  of  free  space 

Nq  and  =  the  vertical  and  horizontal  components,  respec¬ 
tively,  of  the  radiation  vector  in  spherical 
coordinates. 
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Ng  is  multiplied  by 

Kv  =  Kv 3  +  1  -  2  Kv  cos  «rv  -  2  kh  sin  A)  (8.23a) 


and  N  .  by 

0 


XH  =  KHa  +  1  -  2  Kh  cos  (<jfH  -  2  kh  sin  A) 


(8.23b) 


to  give  the  sum  of  the  direct  and  ground  reflected  waves.  Here,  h  is  the 

2tt 

height  above  ground,  k  =  -y  ,  and  X,  and  are  the  vertical  and  horizontal 
ground  reflection  factors,  respectively.  Equations  (8.23a)  and  (8.23b) 
are  valid  for  horizontal  antennas  only. 

8.6.  Antenna  Efficiencies 

The  efficiency  of  an  antenna  or  its  ability  to  radiate  the  power 
supplied  to  it  is  assumed  in  these  methods  to  be  100%  for  horizontal 
non -terminated  antennas,  and  67%  (1.7  dB  reduction)  for  terminated 
horizontal  antennas.  The  efficiency  (F^)  of  an  inverted  "L"  antenna 
(figure  8.  1)  is  represented  by  the  following  polynomial 

Ft  =  20  log  (6.  335  X  +  67.  95  X3  -  693.  00  X3  +  1600.  00  X4) 

10 

X  =  physical  height/  X 

where  (8.24) 

X  =  wavelength  of  the  operating  frequency 

=  theoretical  power  gain  to  be  reduced  by  this  amount-~dB. 

The  efficiency  (Fy)  of  the  grounded  vertical  antenna  (figure  8.  1)  is 
represented  by  the  following  polynomial 

Ft  =  25.  646  -  364.  817  X  +  2179.  89  X2  -  6091. 33  X3  +  6416.  702  X4  (8.25) 
(variables  are  identical  to  above). 


68 


PERCENT  EFFICIENCY 


The  above  efficiencies  are  commonly  used  throughout  the  field 
[Laitinen  1957];  therefore,  no  further  justification  will  be  given  for 
their  acceptance. 

For  a  complete  explanation  of  the  efficiencies  used  and  the  curves 
in  figure  8.  1,  the  reader  is  referred  to  the  work  of  Laitinen  [1957]. 

For  the  reception  of  high  frequency  signals,  the  ability  of  the 
antenna  to  discriminate  between  the  signal  and  the  noise  is  the  important 
consideration,  since  the  major  noise  sources,  atmosphere,  man-made 
and  cosmic,  are  external  to  the  antenna  system.  In  the  absence  of  better 
information,  it  is  assumed  that  this  external  noise  arrives  equally  from 
all  directions,  and  the  receiving  antenna  gain  (performance  of  the  antenna 
in  the  reception  of  this  signal  relative  to  its  performance  in  the  reception 
of  noise)  is  the  directive  gain  of  the  antenna;  i.e.  ,  antenna  efficiency  is 
not  considered  in  the  receiving  antenna  gain. 

The  directive  gain  may  be  applied  to  signal-to-noise  ratio  predictions; 
however,  power  gain  should  be  used  when  predicting  system  loss  or 
receiver  input  powei  . 

9.  Upper  Limits  of  Frequencies 

The  upper  limits  of  fiequencies  are  calculated  from  monthly  median 
values  of  foF2,  the  maximum  height  of  electron  density  which  is  derived 
from  the  M-3000  factors,  and  the  thickness  of  the  parabolic  .layer  (y). 

Conditions  along  the  path  are  estimated  from  the  above  variables 
averaged  from  areas  of  reflection  along  the  path  to  yield  a  typical  profile 
of  height  and  electron  density  associated  with  the  parabolic  layers.  The 
Maximum  Usable  Frequency,  MUF  (0.  50),  is  the  highest  operating  fre¬ 
quency  (f)  that  is  still  reflected  by  the  layer  for  0.  50  of  the  days  within 
the  month  when  the  electron  density,  height,  and  thickness  are  monthly 
median  values.  Only  the  low-angle  ordinary  rays  are  considered.  The 
probability  of  reflection  associated  with  frequencies  equal  to  or  greater 
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than  the  monthly  median  MUF  is  approximated  by  knowing  its  height  of 
reflection,  relation  to  the  MUF  (0.  50)  and  associated  distribution  of  the 
MUF1  s  as  seen  in  the  tables  contained  in  Appendix  B. 

Paths  are  not  considered  to  be  long  (4000  km  or  greater)  or  short 
in  the  calculation  of  the  upper  limits  since  one  hop  may  go  only  3000  km  or 
as  great  as  7000  km  depending  on  operating  frequency,  height  of  maximum 
ionization,  thickness  of  the  layers,  and  the  retardation  (bending)  in  the 
lower  regions  (section  7.  2).  Mixed  layer  MUFs  are  not  included  in  this 
method;  therefore,  slight  discrepancies  may  appear  in  extreme  cases, 
e.  g.  ,  mixed  mode  MUF  (<  55). 

9.  1.  Empirical  Distribution  of  F2(3000)MUF 
An  investigation  of  the  distribution  of  daily  values  of  MUF  about 
their  monthly  median  was  carried  out  in  order  to  provide  information 
needed  to  estimate  "circuit  reliability"  in  the  high  frequency  band.  Three 
points  on  the  distribution  curve  were  considered,  the  values  of  daily  MUF 
exceeded  0.  90,  0.  50,  and  0.  10  of  the  days  within  the  month.  The  ratios 
of  the  decile  values  to  the  median  value,  were  computed  for  comparison 
with  the  often  used  values  of  0.85  and  1.15. 

The  data  used  in  this  study  are  values  of  the  standard  MUF,  the 
product  of  the  foF2  and  M-3000  factor  scaled  from  vertical-incidence 
records.  Data  were  analyzed  from  13  stations  representing  a  range  of 
geomagnetic  latitudes  from  71°S  to  88  °N.  The  temporal  variations  of 
the  MUF-distribution  were  determined  by  considering  observations  at  all 
24  local  time  hours,  each  month  of  the  year,  and  periods  representing 
low,  medium,  and  high  solar  activity. 

The  tables  in  Appendix  B  give  the  ratios  of  upper  and  lower  decile 
MUFs  to  median  MUF.  Each  table  shows  the  values  for  a  given  season, 
a  given  solar  activity,  local  time  hours  00,  04,  08,  etc.  ,  and  each  10° 
of  geographic  latitude  from  10°  to  80°,  north  or  south. 
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The  study  indicated  that  the  distribution  of  MUFs  is  wider  at  night 
than  in  aaytime.  and  wider  at  low  latitudes  than  high  latitudes  in  daytime. 
Again  in  daytime,  the  distribution  is  wider  in  summer  than  winter,  except 
at  high  latitudes  where  the  reverse  is  true.  The  sunspot  number  depend¬ 
ence  is  weaker,  but  in  daytime  the  difference  between  the  two  ratios  seems 
to  increase  with  sunspot  number  at  latitudes  higher  than  40°  and  to  de¬ 
crease  with  increasing  sunspot  number  at  latitudes  below  40°  [Davis  and 
Groome  1965].  The  study  indicated  that  the  distribution  was  mostly  a 
function  of  the  foFZ  and  not  the  M(3000)F2;  therefore,  the  distributions 
are  assumed  valid  for  any  oblique  path. 

9.2.  Empirical  Distribution  of  Transmission  Loss  (Yp) 

The  day-to-day  variations  of  signals  on  operating  circuits  of 
various  lengths--55  to  15,  689  km  in  the  temperate  region  (below  40° 
geomagnetic  latitude)  and  1220  to  4470  km  in  the  vicinity  of  the  auroral 
zone  (above  50°  geomagnetic  latitude) --have  been  analyzed.  The  paths  were 
inspected  at  a  given  hour  to  best  determine  those  signals  being  propagated 
by  frequencies  below  the  daily  MUF.  The  results  showed  that  the  distrib¬ 
utions  were  distinct  functions  of  geomagnetic  latitude,  season,  local  time 
and  length  of  path.  The  distributions  of  transmission  loss  are  expressed 
in  terms  of  the  differ 'nee  from  the  median  of  the  lower  decile  (D  )  and 

i/ 

the  upper  decile  (Du).  This  was  necessary  because  the  distributions  were 
not  log-normal  but  were  assumed  to  be  log-normal  on  either  side  of  the 
median  values.  The  tables  in  Appendix  C  were  condensed  from  the  work 
of  Davis  and  Groome  [1964],  Laitinen  and  Haydon  [1962]  and  Haydon  and 
Lucas  [1966],  They  show  the  separation  from  the  median  of  the  trans¬ 
mission  loss  exceeded  84  percent  of  the  time  (S^  =  D^/  1.  28)  and  the 
transmission  loss  exceeded  16  percent  of  the  time  (Su  =  Du  /  1.28). 

A  study  of  the  distributions  shows  that  the  losses  tend  to  be  greater 
for  paths  in  the  region  of  64°  geomagnetic  latitude  and  more  widely 
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distributed  with  the  diurnal  maximum  falling  between  04-10  local  mean 
time.  The  seasonal  maximum  loss  occurs  at  equinox  for  both  long  and 
short  paths,  with  the  short  path  generally  showing  greater  excess 
attenuation  at  high  latitudes.  The  distributions  are  probably  less  reliable 
on  paths  near  the  geomagnetic  equator  since  most  of  the  temperate  paths 
which  were  analyzed  fell  above  15°  north  geomagnetic  latitude. 

10.  Reliability,  q(Rh) 

This  section  will  be  mainly  concerned  with  suggesting  a  method  of 
estimating  the  probability  of  an  adequate  signal-to-noise  ratio  at  a  given 
hour,  i.e.  ,  the  circuit  reliability  of  HF  sky-wave  paths.  "Circuit 
Reliability"  is  defined  as  the  probability  of  successful  communications 
at  a  given  hour  within  the  month  at  a  specific  operating  frequency. 

A  consideration  in  the  evaluation  of  a  sky-wave  circuit  is  the 
maximum  frequency  that  will  be  supported  at  a  given  hour  or  conversely 
what  is  the  probability  that  a  chosen  frequency  will  be  propagated. 

Tables  of  Appendix  B  show  the  distribution  of  the  F2  region  Maximum 
Usable  Frequencies  as  a  function  of  time  of  day,  season,  solar  activity, 
and  geomagnetic  latitude  of  the  circuit.  These  distributions  of  the 
F2  layer  Maximum  Usable  Frequencies  (MUF)  are  daily  variations  about 
the  monthly  median  MUF  as  predicted  by  the  Institute  for  Telecommuni¬ 
cation  Sciences  and  Aeronomy  (formerly  CRPL). 

The  distribution  of  the  MU?'  on  either  side  of  the  median  is  assumed 
to  be  log-normal;  thus  the  probability  of  support  at  a  chosen  frequency 
propagated  at  a  given  elevation  angle  can  be  determined  by 


q# 
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(10.1) 
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where  f  =  operating  frequency  -  MHz 
aMU r  ^  =  (9th  or  1st  decile)/  1.28 

wuf  £  =  frequency  supported  0.  50  of  tne  time  at  angle  A  . 

The  above  is  only  an  approximation  since  the  distribution  is  skewed 
as  can  be  seen  in  the  tables  of  Appendix  B.  In  practice,  however,  the 
errors  associated  with  the  log-normal  assumption  (either  side  of  vhe 
median)  may  be  considered  to  be  negligible  in  typical  communication 
problems. 

The  frequency  which  has  a  0.  9  or  a  0.  1  probability  of  sky-wave 
propagation  at  a  given  hour  within  the  month  can  be  estimated  by  using 
the  appropriate  decile  values  with  the  predicted  monthly  median  MUF 
provided  the  lower  region  retardation  and  depth  of  penetration  into  the 
reflecting  region  does  not  significantly  change  the  takeoff  angle  (A). 

The  frequencies  having  a  0.  9  and  a  0.  1  probability  of  sky-wave 
propagation  correspond  to  the  classically  defined  Optimum  Traffic  Fre¬ 
quency  (FOT)  and  the  Highest  Probable  Frequency  (HPF). 

The  probability  of  a  given  frequency  being  propagated  by  a  given 
sky-wave  mode  can  be  estimated  from  the  distributions  of  the  critical 
frequencies,  layer  heights,  and  semithicknesses  of  the  ionospheric 
regions  associated  with  the  sky-wave  mode. 

Successful  communications  also  depend  upon  an  adequate  signal-to- 
noise  ratio  at  the  receiver.  Variations  of  the  ionospheric  and  geophysical 
parameters  result  in  a  variation  of  the  received  signal.  This  variation 
can  be  expressed  as  a  distribution  about  the  monthly  median.  Tables  of 
Append'-  c  C  show  this  distribution  as  a  function  of  the  season,  the  time  of 
day,  the  geomagnetic  latitude  and  the  length  of  the  circuit  for  frequencies 
below  the  daily  MUF.  The  probability  of  successful  HF  sky-wave  commu¬ 
nications  depends  upon  the  probability  of  the  operating  frequency  being 
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supported  via  the  ionosphere  and  the  probability  that  the  signal-to-noise 
ratio  will  exceed  some  acceptable  level. 

These  events  are  assumed  to  be  independent;  therefore,  the  com¬ 
bined  probability  is  the  product  of  the  probability  that  the  signal-to-noise 
ratio  exceeds  a  given  level  and  the  probability  that  at  least  one  mode  will 
be  present  to  produce  the  median  available  signal-to-noise.  The  total 
probability  of  ionospheric  support  at  a  given  frequency  f  is  taken  to  be 
the  ray  path  producing  the  highest  probability  at  a  given  hour.  The 
probability  q^  is  evaluated  by  the  following  integral  assuming  no 
correlation  between  the  received  signal  and  the  received  noise  powers. 
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(10.2) 


where 

Rh  =  minimum  acceptable  signal-to-noise  ratio  -  dB 
s  =  monthly  median  signal  -  dB 

TT  =  monthly  median  of  the  hourly  median  radio  noise  -  dB 
aH  =  standard  deviation  of  radio  noise  -  dB 
as  =  standard  deviation  of  received  signal  -  dB. 


The  resulting  probability  of  exceeding  a  minimum  required  signal- 
to-noise  ratio  on  a  given  day  is  then  q(Rh)  =  (q*)  •  (qs/N).  This  total 
probability  is  defined  as  the  fraction  of  the  days  within  the  month  at  the 
given  hour  that  the  requirement  (Rh)  will  be  met  (CIRCUIT  RELIABILITY). 

The  median  noise  power  (N)  and  its  associated  distribxitions  (<JN)  are 
calculated  using  numerical  representations  of  CCIR  322  [Lucas  and  Harper 
1965],  The  median  excess  signal  power  (5)  above  the  quasi-maximum 
signal  power  and  its  distribution  (CJS)  are  found  in  Appendix  C  [Davis  and 
Groome  1965]  and  [Laitinen  and  Hay  don  1962]. 
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The  above  method  can  be  applied  to  propagation  via  the  regular  E, 
sporadic  E  or  FI  layers  as  well  as  the  F2  layer  if  sufficient  knowledge 
of  the  distributions  is  available  and  independence  can  be  assumed. 

The  variation  of  the  regular  E  layer  MUF'  s  is  very  small  and 
therfore  a  qf  of  1  can  be  assumed  when  f  is  less  than  the  E-MUF  and 
a  qf  of  0  when  f  is  greater  than  E-MUF.  The  signal  distributions  are 
assumed  to  be  valid  for  both  E  layer  and  F2  layer  propagation  since  in  the 
empirically  determined  distributions  no  attempt  was  made  at  sorting  out  the 
individual  ray  paths  of  propagation.  The  contribution  of  each  individual 
ray  path  of  propagation  is  in  part  contained  in  the  distributions  about  the 
median  values. 

It  is  possible  using  the  above  approach  for  the  engineer  or  circuit 
planner  to  inspect  the  results  of  the  above  equations  and  extract  much 
more  useful  information  than  in  the  past.  For  instance,  one  could  look 
at  qsyN  and  determine  if  a  change  in  power,  antenna  or  modulation  would 
yield  a  better  overall  circuit  reliability  due  to  an  increase  in  available 
signal-to-noise  or  a  decrease  in  the  basic  required  signal-to-noise  ratio. 
The  quantity  qf  will  dictate  the  frequency  at  which  an  increase  of  power, 
etc.  ,  would  be  advisable,  since  values  of  qf  dictate  the  overall  circuit 
reliability  at  frequencies  which  yield  low  values  of  qf. 

Conversely,  high  values  of  qf  at  a  given  frequency  dictate  that  an 
increase  of  available  signal-to-noise  ratios  or  a  decrease  in  required 
signal-to-noise  may  yield  an  increase  ir.  the  overall  ciicuit  reliability. 

A  careful  inspection  of  tables  10.  1,  10.2  and  10.  3  will  show  some 
possible  effects  of  power  and/  or  antenna  changes  at  given  frequencies  on 
the  overall  "circuit  reliability"  of  a  typical  high  frequency  communication 
circuit.  (See  section  15  for  definition  of  variables  appearing  in  the  above 
tables. ) 

For  example,  at  0200  GMT  using  an  operating  frequency  of  9  MHz 
a  power  increase  to  30.00  kW  (Table  10,  1}  from  1.00  kW  (Table  10.2) 
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Table  10.  1. 
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increases  the  circuit  reliability  from  .  72  to  .  98.  An  inspection  of  the 
same  frequency  at  the  same  hour  in  Table  10.3  indicates  that  a  change 
of  antennas  from  curtains  and  rhombic s  to  a  horizontal  dipole  and  an 
assumed  of  3  dB  lowers  the  circuit  reliability  to  .  07  when  using 
l  kW  of  pow:  r. 

The  prediction  errors  involved  in  the  estimation  of  the  radio 
noise  levels  and  the  available  signal  levels  have  been  neglected  in  the 
calculation  of  the  circuit  reliability.  It  is,  therefore,  evident  that  any 
value  of  circuit  reliability  will  under  these  circumstances  be  the  most 
likely  value  of  the  parameter  with  actual  values  lying  on  either  side  with 
a  like  probability  („  5)  if  the  frequency  has  ionospheric  support.  The 
prediction  errors  associated  with  the  expected  time  availability  qP  are 
included  m  the  following  sections  dealing  with  service  probability  (Qr). 

The  concept  of  service  probability  was  originally  developed  by  K.  A.  Norton 
in  1958  in  connection  with  predicting  the  performance  of  tropospheric 
communication  circuits. 


11.  Grade  of  Service  (gr)  and  Required  Signal-To-Noise  Ratios  (Rh) 

The  "Grade  of  Service  (gr  )  refers  to  the  degree  of  reliability  over 
a  short  period  of  time  (approximately  a  few  minutes  to  an  hour)  during 
which  the  statistics  of  the  signal-to-noise  ratio  may  be  considered  to  be 
stationary.  It  can  be  expressed  for  example  as  the  percentage  of  error- 
free  messages,  the  intelligibility  achieved,  or  the  percentage  of  satisfied 
observers.  "  [CCIR  Report  322] 

The  required  signal-to-noise  ratio  (Rh)  at  the  receiver  for  a  given 
grade  of  service  (gr )  will  have  some  variation  due  to  variations  of  the 
short-term  noise  and  signal,  multipath  propagation,  circuit  equipments, 
etc.  There  is  a  prediction  error  ar  in  defining  the  required  signal-to- 
noise  ratio  for  any  given  grade  of  service.  This  report  assumes  that  this 
prediction  error  is  two  decibels  for  all  grades  of  service  for  the  lack 
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of  better  information.  If  the  circuit  engineer  has  information  on  the 
prediction  errors  involved  in  the  estimate  of  the  required  signal-to-noise , 
it  should  be  used  in  lieu  of  the  suggested  two  decibels. 

12.  Time  Availability  (qr )  and  Service  Probability  (QT) 

"The  time  availability  refers  to  the  percentage  of  the  hours  or 
other  short  periods  of  time  used  in  defining  the  grade  of  service,  during 
which  the  specified  grade  of  service  or  better  is  achieved.  The  time 
involved  should  include  all  of  the  expected  variations  and  may  be  an 
entire  sunspot  cycle,  a  year,  a  particular  season  or  month,  or  certain 
hours  of  the  day  during  a  specified  longer  period.  "  [CCIR  Report  322] 
This  specified  longer  period  is  the  given  hours  of  the  day  within  a  given 
month  of  the  year.  This  means  that  when  0.90  time  availability  is  speci¬ 
fied  a  given  grade  of  service  (required  signal-to-noise  ratio)  is  expected 
to  be  exceeded  0.  90  of  the  days  within  the  month  for  the  specified  hours. 

"Service  Probability  (Qt  )  is  defined  as  the  probability  that  the 
specified  grade  of  service  (gr)  or  better  will  be  achieved  for  the  specified 
time  availability  (qr)."  (CCIR  Report  322)  This  combines  statistically  the 
uncertainties  of  the  many  parameters  involved  in  the  calculation  of  the 
performance  of  high  frequency  circuits. 

The  uncertainty  (aR  )  of  predicting  a  required  signal-to-noise  ratio 
for  a  specified  grade  of  service  has  been  explained  in  section  11.  The 
uncertainties  (op  ,  o3u  and  a3<^)  of  predicting  the  median,  upper  decile  and 
lower  decile  respectively  of  system  loss  or  conversely  the  available 
signal  level  at  the  receiver  are  shown  in  tables  1  through  3  of  Appendix  D. 
The  values  of  prediction  errors  included  in  Appendix  D  were  obtained 
from  comparisons  of  predicted  values  of  received  signal  with  those 
observed  on  operating  circuits.  The  comparisons  include  analyzed 
circuits  ranging  in  great-circle  distance  from  250  km  to  16,000  km  and 
using  frequencies  ranging  from  2.  06  MHz  to  20.  0  MHz.  A  wide  range 
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cf  latitudes  was  included  in  the  circuits  analyzed  with  most  circuits 
lying  between  12  and  80  degrees  north  geomagnetic  latitude.  Sufficient 
data  were  not  available  for  high  and  low  solar  activity  levels  to  attempt 
a  correlation  with  solar  activity.  An  attempt  was  made  to  include  a 
frequency  dependence  in  the  uncertainties;  however,  no  consistent 
variation  was  apparent.  A  distinct  variation  did  appear  as  a  function 
of  local  time  of  day  and  geomagnetic  latitude  of  the  circuit.  The  values 
of  Appendix  D  are  a  function  of  these  two  variables  and  are  assumed  to 
be  normally  distributed  either  side  of  the  median  ap. 

The  uncertainties  of  predicting  the  radio  noise  level  at  the  receiving 
site  are  contained  in  CCIR  Report  322  as  the  error  in  predicting  (1)  the 
median  1  MHz  atmospheric  noise  level  ) ,  (2)  the  lower  decile  value 

of  the  atmospheric  noise  l-;vel  ),  and  (3)  the  upper  decile  value  of 
the  atmospheric  noise  levs!  (o0u).  The  three  variables  appearing  above 
are  evaluated  using  the  numerical  coefficients  of  tables  1  through  3  of 
Appendix  E  in  the  following  equation 

Y(x)  =  Ax  +  AsX  +  A3X2  +  A*X3  +  AgX4  {12.  1) 

where 

Y(x)  =  arAM,  a0L  or  q,u 

X  =  loge  (f),  and 

f  =  operating  frequency  -  MHz 

12.  1  Calculation  of  the  Protection  Factor  (C)  for  a  Specified 

Time  Availability  (qr  ) 

The  day-to-day  variations  of  the  hourly  median  signal  levels  (S) 
and  the  hourly  median  noise  levels  (N)  are  assumed  to  be  normally  dis¬ 
tributed  either  side  of  the  median  values  and  can  be  described  by  the 
median  values  and  the  deviation  Ds  N,  of  the  value  exceeded  some  speci¬ 
fied  fraction  of  the  time,  from  the  median.  A  protection  factor  (C) 


needed  to  provide  the  required  signal-to-noise  ratio  (Rh)  for  the  specified 
time  availability  (qr )  can  be  determined- -assuming  no  correlation- -from 


C2  =  Dn2  +  Ds2  (12.2) 

Dn  =  values  of  the  average  noise  power  exceeded  with 
probability  (1  -  qp) 

Ds  =  value  of  the  average  signal  power  exceeded  with 
probability  qr. 


The  values  DN  and  Ds  are  evaluated  using  the  values  of  tables  1 
through  12  of  Appendix  C  and  the  numerically  represented  values  of 
FAH  ,  DL  ,  and  D„  of  CCIR  Report  322  by  the  evaluation  of  the  function 


(12.3) 


then 


a.  =  given  time  availability  (0  ^  qr  ^  1) 

A  =  standard  normal  deviate  corresponding  to  a  specified 
value  of  qr 


Ds  =  A  Su  and 
Dn  =  a  Du 


(12.4) 


The  values  of  time  availability  are  assumed  to  be  between  0.  50  and 
1  in  the  above  illustration;  however,  small  changes  in  the  above  equations 
will  allow  the  calculations  to  be  performed  for  values  of  time  availability 
less  than  0.  50. 


12.  2.  Calculation  of  Service  Probability  (Or  ) 

The  service  probability  (QT  )  corresponding  to  the  specified  time 
availability  (qr)  for  a  given  required  signal-to-noise  ratio  (Rh  )  can  be 
calculated  using  the  value  of  C  calculated  previously  and  the  required 
sienal-to-noise  ratio  (Ru)  of  section  11  and  the  total  uncertainties  (a, . )  to 
be  explained  later. 
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The  median  value  of  the  expected  required  signal  level  S(qr)  is 
evaluated  by 

S(qP)  =  Fam  +  C  +  Rh  (dB  <  1  watt)  (12.  5) 

i'AM  =  median  amplitude  of  the  hourly  median 
of  the  seasonal  median  prevailing  radio 
noise  level 

C  =  required  protection  factor  for  the  specified 
time  availability  (qP). 

Rn  =  required  signal -to-noise  ratio  for  the 
specified  grade  of  service  (gr). 


The  total  uncertainties  (ct)  involved  in  the  prediction  of  circuit 
performance  can  be  evaluated,  on  the  assumption  the  prediction  errors 
are  uncorrelated,  from 

o,2  =  q2  +  ap2  (s)  +  ac2  (n)  +  a„2  (s)  +  (dB)2  !12’6) 

ap  =  standard  error  in  decibels  in  predicting  the  required 
signal-to-noise  ratio  (Rh  ) 

Op  =  root  mean  square  error  in  decibels  in  predicting 
the  median  signal  level  (S) 

oFAM  =  root  mean  square  error  in  decibels  in  predicting 
the  median  radio  noise  level  (N) 

q.  (n)  =  additional  standard  error  in  decibels  m  predicting 
the  radio  noise  level  for  a  given  median  level,  at 
the  given  value  of  time  availability  (qP) 

oe(s)  =  additional  standard  error  in  decibels  in  predicting 

the  signal  level  for  a  given  median  level,  at  the  given 
value  of  time  availability  (qr). 


Since  the  value  of  the  expected  required  signal  level  S(qP  )  has  been 
evaluated  along  with  the  total  uncertainties  of  the  predictions  (Ot),  the 
service  probability  can  be  evaluated  for  those  frequencies  having  a 
probability  of  support  (q*  )  equal  to  1  by 
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Qf  =  service  probability  for  frequencies  sure 
of  support  via  the  ionosphere 

S  =  predicted  available  median  signal  level- 
dB  <  1  watt 

S(qr)  =  the  median  value  of  the  expected  required 
signal  power  (for  the  given  time  avail - 
ability-dB  <  1  watt) 

crt  =  total  prediction  uncertainties  of  the  signal 
and  noise  levels -dB. 


For  those  frequencies  having  a  probability  of  support  (qf ) 
less  than  one  the  total  service  probability  would  be 

Or  =  Or*-  qf  (1Z„ 


qf  =  probability  of  ionospheric  support 
of  the  operating  frequency  (f) 

Qr  =  service  probability  including 

the  like  lihood  of  support  via  the 
ionosphere. 

During  the  hours  at  a  given  operating  frequency  when  the  atmos¬ 
pheric  radio  noise  is  not  prevailing,  the  following  values  of  the  predic¬ 
tion  errors  for  cosmic  or  man-made  noxse  are  recommended  if  more 
reliable  information  is  not  empirically  available. 


Cosmic  Noise  (dB)  Man-Made  Noise  (dB) 

Qf  AM  °0u  ^OL  QfK*  <%»  °DL 

0.  5  0.  2  0.  2  3.  0  1.  5  1.5 

Tables  12.  1  and  12.  2  show  predictions  of  service  probability  for 
a  chosen  time  availability  of  0.  90  and  0.  80.  respectively. 
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Sample  Computer  Print-Out  Showing  Service  Probability 
for  a  Chosen  Time  Availability  of  0,  90  Using 
Described  Prediction  Model 
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for  a  Chosen  Time  Availability  of  0.  80  Using  Described 

Prediction  Model 
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An  inspection  of  the  above  tables  will  show  that  the  service  prob¬ 
ability  (Or)  will  be  dominated  by  the  probability  of  support  (qf )  at  the 
higher  frequencies;  but  will  be  controlled  by  the  available  signal-to-noise 
ratio  and  the  prediction  errors  (Ot )  at  the  lower  frequencies. 

13.  Multipath  Considerations 

The  usual  type  of  multipath  which  provides  the  greatest  difference 
in  propagation  times  results  from  two  or  more  paths  with  a  different 
number  of  hops.  Multipath  propagation  such  as  this  can  occur  a  large 
percent  of  the  time  over  a  wide  frequency  range.  It  can,  however,  be 
minimized  by  operating  at  frequencies  where  the  likelihood  of  multiple 
signals  of  comparable  strength  is  small,  e.g.  ,  as  close  to  the  MUF  as 
practical. 

It  is  assumed  that  an  adequate  signal-to-noise  ratio  will  suffice  to 
determine  circuit  reliability  if  the  operator  has  an  adequate  choice  of 
frequencies  in  the  complement  to  permit  operation  where  the  likelihood 
of  multiple  signals  is  low'.  The  multipath  predictions  are  intended  to  be 
used  as  an  aid  to  determine  the  likelihood  of  multipath  at  a  given  frequency 
at  a  specific  hour.  If  the  reliability  or  service  probability  of  two  fre¬ 
quencies  is  equal,  but  one  possesses  a  greater  likelihood  of  multipath, 
this  will  also  assist  in  the  choice  of  operating  frequency. 

The  difference  in  propagation  times  for  various  ray  paths  is  a 
function  of  frequency  end  transmission  distance;  therefore,  the  probability 
of  multipath  is  also  dependent  upon  these  variables.  The  multipath  which 
is  predicted  is  based  upon  some  discrete  minimum  difference  in  received 
power  and  maximum  tolerable  time  delay  between  ray  paths.  The  prob¬ 
ability  of  multipath  is  assumed  to  be  the  probability  that  a  secondary  ray 
path  will  be  present  to  produce  minimum  difference  in  received  signal 
powers. 
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In  addition  to  selecting  the  frequency  which  is  less  likely  to  be 
subject  to  multipath,  the  multipath  predictions  may  also  be  applicable  to 
antenna  design  or  selection  to  minimize  multipath. 

It  is  likely  that  predictions  of  this  type  will  become  more  important 
as  higher  data  rate  systems  are  employed.  The  predictions  are  designed 
to  estimate  multipath  probability  as  a  function  of  time  delay  and  minimum 
permissible  level  of  the  delayed  signal.  Tables  15.2a  and  15.2b  show 
some  typical  predictions  of  frequencies  which  possess  multipath.  They 
have  an  asterisk  above  the  mode  designation  and  an  associated  multipath 
probability  printed  below  the  circuit  reliability  or  service  probability. 

14.  Comparisons  of  Predictions  with  Observations 

Limited  comparisons  of  the  predictions  with  observations  are 
shown  to  illustrate  ihe  type  of  results  which  can  be  expected  (figures  14.  1 
through  14.  8). 

Absolute  comparisons  between  predictions  and  observations  are 
normally  difficult  because  the  observed  field  strengths  are  usually 
abbreviated,  e.g.  ,  read  on  a  receiver  ”S"  meter.  Also  changes  in 
system  sensitivity  with  frequency  and  insufficient  knowledge  of  the  actual 
equipment  (antennas,  transmitter  power,  line  losses,  multicouplers,  etc.) 
are  common. 

The  checks  shown  have  some  of  the  difficulties  mentioned  above, 
and  assumptions  were  required  in  many  of  the  predictions.  When  the 
details  of  the  equipment  parameters  and  methods  of  recording  were 
not  sufficiently  known,  only  relative  values  are  shown. 

Figure  14.  1  shows  the  agreement  of  observed  monthly  median 
maximum  observed  frequency  MOF  (0.  50)  with  the  predicted  monthly 
median  MUF  (0.  50).  While  the  curves  depart  somewhat  at  certain  hours 
the  diurnal  shape  is  in  quite  good  agreement.  The  root  mean  square  error 
of  the  predicted  curve  in  figure  14.  1  was  calculated  to  be  0.  94  MHz. 


B9 


GREENWICH  MEAN  TIME 


Figure  14.  1.  Fort  Monmouth  -  Palo  Alto  Smooth  MOF  (0.  50)  Data 
Versus  Predicted  MUF  (0.  50)  (November  SSN  65) 
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Figure  14„2.  Predictions  and  Observations  of  the  Elmendorf 
McClellan  Oblique  Sounder  January  1963  -  12  kW  -  SSN  29 
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Figure  14.  5,  WWV  Monitoring  V*  rsus  Predictions  Beltsville 
Crowborough  (June  -  SSN  14-10  MHz) 
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Figure  14.8.  Predicted  Versus  Observed  Relative  Received 
Voltage  Taipei  -  Tokyo  -  1145  N.  M.  - 
June  -  SSN  161  -  200  Watts  -  Rhombics 
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Figure  14.2  shows  a  comparison  of  predicted  MUF  (0.  5)  and 
oblique  sounder  MOF  (0.5)  [Williams  and  Egan  1963].  The  comparison 
may  not  be  optimum  since  the  predicted  MUF  should  agree  more  closely 
with  the  junction  frequency  as  recorded  on  the  oblique  sounder.  The 
comparison  was  made,  however,  since  both  quantities  are  a  measure 
of  the  upper  limit  of  fiequencv.  The  root  mean  square  error  of  the 
predicted  MUF  (0.  5)  was  calculated  to  be  1.  55  MHz.  The  root  mean 
squared  error  of  the  predicted  LUF  (0.50)  is  .43  MHz. 

Figures  14.3  through  14.8  are  comparisons  of  observed  received 
signal  power  with  predicted  received  signal  power.  The  figures  shown  are 
representative  of  many  years  of  data  analyzed  at  ITSA  (formerly  CRPL) 
and  standard  errors  of  the  estimates  have  not  been  calculated  for  the 
specific  figures  shown.  The  agreement  between  observations  and  pre¬ 
dictions  appears  good  except  for  possibly  figure  14.  5,  where  the  predicted 
received  power  for  the  middle  of  the  day  is  in  poor  agreement  with  the 
observations.  No  really  satisfactory  explanation  for  this  is  currently 
available  except  that  it  is  a  high  latitude  path  with  a  high  probability  of 
sporadic  E  and  current  prediction  methods  do  not  include  sporadic  F. 

15.  High  Frequency  Sky-Wave  Computer  Prediction  Routine 

The  computer  program  is  written  primarily  in  Fortran  3600  for 
the  CDC-3600  machine.  The  program  is  self-contained  except  for  a  few 
library  subroutines  which  are  assumed  to  be  available  on  most  machines. 
Comment  cards  are  included  throughout  the  program  to  provide  easy 
stepping  stones  for  those  familiar  with  computer  programming.  All 
statements  are  intended  to  be  simple  and  straightforward,  utilizing  no 
peculiarities  of  a  specific  monitor  or  machine .  The  basic  Fortran 
should  allow  simplified  conversion  to  other  machines  and  languages  if 
desired.  All  parameters  which  are  evaluated  many  times  are  included 
as  subprograms.  The  program  is  fast;  yet  ease  of  incorporating 
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changes  and  new  knowledge  is  not  sacrificed  for  speed  alone.  The 
program  is  intended  to  be  a  basic  program  on  which  to  build  as  new 
knowledge  and  methods  become  available  from  basic  research. 

Tables  15.1  through  15.8  and  figures  15.1  and  15.2  show  the  types 
of  output  possible  using  the  basic  concepts  explained  in  preceding  sections. 
Many  items  on  the  printed  sheets  are  self-explanatory,  but  for  the  sake 
of  completeness,  the  explanation  is  as  follows: 

A.  Heading  reading  left  to  right 

First  Line 

1.  sequential  number  of  circuit  as  entered  into  machine  (X); 

2.  month  of  year  (AAA); 

3.  solar  activity  level  (SSN  =  XXX); 

4.  customer  identification  number  derived  from  the  great- 
circle  distance  and  backward  bearing  (AA  XX.  XXX). 

Second  Line 

1.  transmitter  and  receiver  locations; 

2.  azimuths  (forw'ard  and  backward); 

3.  nautical  miles  (N.  Miles). 

Third  Line 

1.  geographic  locations  of  transmitter  and  receiver  in 
hundredths  of  degrees  (XX.  XXA-XXX.  XXA); 

2.  azimuths  (XXX.  X  XXX.  X); 

3.  great-circle  distance  of  path  in  nautical  miles  (XXXX.  X). 

Fourth  Line 

1.  transmitting  antenna  and  description; 

2.  receiving  antenna  and  description. 

Note:  A  minus  number  preceding  any  descriptive  symbol  denotes 
aparameter  in  wavelengths,  i.e.  ,  -12  =\/Z,  -2  =  2X., 

-14  =  X/4.  All  other  symbols  are  self-explanatory. 
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Fifth  lane 

1.  off  azimuth  of  transmitting  antenna  (XXX  degrees); 

2.  minimum  angle  above  the  horizon  for  which  any  mode 
will  be  calculated  (XX  degrees); 

3.  off  azimuth  of  receiving  antenna  (XXX  degrees). 

Sixth  Line  (if  any) 

1.  output  power  of  the  transmitter  at  the  antenna  terminals; 

2.  (some  options) 

measured  or  assumed  man-made  noise  level  at  the 
receiving  antenna  site  (d'B  <  1  watt  at  3  MHz  in  a  1  Hz 
bandwidth); 

3.  (some  options) 

the  Lowest  Useful  Frequency  (minimum  operating  fre¬ 
quency),  i.  e.  ,  0.  50  indicates  lowest  frequency  having 
a  circuit  reliability  of  0.  50 ; 

4.  required  signal-to-noise  ratio  (dB) 

a.  the  median  signal  power  required  in  the  occupied 
bandwidth  relative  to  the  noise  power  in  a  1  Hz 
bandwidth ; 

b.  the  signal  power  required  must  be  the  same  type 
of  power  indicated  at  the  transmitting  antenna 
terminals. 

Seventh  Line  (if  any) 

1.  minimum  tolerable  difference  in  amplitude  of  received 
signal  via  two  different  ray  paths,  (M.P.PWR) 

2.  maximum  tolerable  delay  time  of  two  different  ray  paths 
(M.  P.  DELAY). 

Note;  The  criteria  of  1  and  2  above  must  both  be  met  for 
the  same  two  modes  for  detrimental  multipath  to  be 
indicated. 
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Eighth  Line  (if  any) 

1.  Greenwich  mean  time  -  hours; 

2.  monthly  median  Maximum  Usable  Frequency  -  MHz 
(XX.  X); 

3.  operating  frequency  complement  -  MHz; 

4.  Time  avail;  specified  time  availability  (fraction  of  days). 

B.  Body  of  print-outs 

MODE:  the  sky-wave  path  contributing  most  to  the  overall 
probability  or  fraction  of  days  a  sky-wave  path  is 
expected,  i.  e.  ,  IE  =  one  hop  via  regular  E,  2X  =  one 
hop  via  F2  and  one  hop  via  regular  E; 

ANGLE;  the  take  off  and  arrival  angle  associated  with  the 
most  probable  sky-wave  path; 

F.DAYS;  fraction  of  the  days  that  at  least  one  sky-wave 
mode  is  likely  to  be  present  to  produce  the 
quasi -minimum  loss; 

LOSS.  DB;  system  loss  as  defined  in  section  5.3; 

DBU;  the  incident  field  intensity  at  the  receiving  antenna 
site,  dB  relative  to  1  uv/ meter; 

DBW;  the  power  at  the  receiving  antenna  terminals,  dB 
relative  to  1  watt; 

S/N;  tie  signal-to -noise  ratio  at  the  receiving  antenna 
terminals  -  dB; 

REL:  circuit  reliability  as  explained  in  section  10; 

5.  PROB;  service  probability  as  explained  in  section  12; 

FOT;  Optimum  Traffic  Frequency  based  on  empirical 

distributions  of  the  MUF  as  explained  in  section  10; 

DELAY:  transmission  time  of  most  probable  modes, 
tenths  of  milliseconds; 

F.  DAYS  M.  P. :  probability  of  multipath  or  fraction  of  days 
multipath  is  expected  (see  section  13). 
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Additional  Comments 


1.  A,  -3,  appearing  in  any  LUF  column  or  row  denotes 
a  LUF  below  3  MHz. 

2.  An  X  appearing  on  the  bottom  line  of  any  graphical 

print -out  indicates  a  reliability,  equal  to  or  exceeding  the 
required  LUF  level,  did  not  exist  at  any  frequency 
{3-30  MHz). 

3.  A  column  below  any  frequency  containing  all  minus  signs 
indicates  that  the  probability  of  all  sky-wave  paths 
inspected  fell  below  0.  05. 

The  computer  program  requires  only  one  data  tape  for  all  needed 
data  and  one  auxiliary  tape  for  "scratch"  purposes.  Appendix  F  containing 
the  program  and  instructions  is  available  from  ITSA.  The  following  items 
are  either  included  or  available: 

1.  Data  tape  written  in  high  (556  ch/  in)  or  low  (200  ch/  in) 
density  in  either  the  binary  coded  decimal  or  binary  format. 

2.  Hollerith  deck  of  complete  program  including  all  subprograms 
and  machine  binary  library  routines  if  necessary. 

3.  Compiled  listing  and  sample  circuits  run  from  compiled 
deck. 

4.  Complete  explanation  of  all  input  and  output  r  uuire<.sents. 

5.  Machine  listing  and  compiled  program  if  necessary. 

6.  Pictorial  views  of  input  decks  and  output. 

7.  Arrangement  of  coefficients  on  data  tape. 

8.  Complete  listing  of  data  tape. 

9.  The  appendix  will  always  contain  latest  improvements  with 
complete  explanation  of  same. 

The  program  can  be  obtained  by  writing: 

Frequency  Utilization  Section 
Ionospheric  Telecommunications  Laboratory 
Institute  for  Telecommunication  Sciences  and 
Aeronomy 

ESSA./NBS  Laboratories 
Boulder,  Colorado  30302 
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Note:  Updating  of  input  va  es  is  constantly  taking  place;  therefore,  values 
in  this  table  and  tab.es  to  follow  are  only  tentative.  Absolute  values  to 
be  checked  will  always  appear  in  Appendix  F. 


Sample  Computer  Print-Out  of  MUF  (0,  50)  and  FOT 
Only  Using  Described  Prediction  Model 
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Figure  15.  1.  Sample  Graphical  Computer  Print-Out  of  MUF  (0.  50) 
and  FOT  Using  Described  Prediction  Model 
(London-Madrid.  December  -  SSN  75) 
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Fraction  of  Days  Multipath  Using  Described  Prediction  Model 
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16.  Conclusions 


16.  1.  Long  Term  Predictions 

World-wide  prediction  of  the  necessary  ionospheric  and  geophysical 
parameters  icnot  considered  to  be  available  currentlv  in  sufficient  detail 
to  merit  sophisticated  ray  tracing  solution  to  typical  problems  in  the 
long  term  prediction  of  the  performance  of  ionospheric  telecommunications 
circuits.  The  philosophy  of  basing  high  frequency  predictions  on  the 
equivalence  theorem,  secant  law,  and  typical  electron  density  profiles 
generated  from  available  ionospheric  predictions  is  generally  superior 
to  simple  ray  path  geometry  between  the  earth  and  concentric  ionospheric 
shells  using  the  classical  control  point  methods. 

The  empirical  relationships  included  in  this  report  should  be  con¬ 
sidered  as  interim  pending  the  development  of  more  refined  relationships 
by  basic  researchers.  It  is,  therefore,  apparent  that  more  confidence 
should  be  placed  in  relative  values  produced  by  prediction  methods  such 
as  the  one  described  than  m  the  absolute  magnitudes  of  these  predictions. 

Caution  should  always  be  taken  when  using  values  produced  by  any 
prediction  routine  based  on  morphology  of  prediction  parameters.  The 
predictions  are  considered  to  be  especially  useful  for  long  term  plannin 
such  as  circuit  design,  siting,  antenna  considerations,  and  frequenc 
assignments.  Additional  work  on  the  correlation  of  past  circuit  per¬ 
formance  with  "numerically"  mapped  parameters  scaled  from  vertical 
ionograms  on  a  world-wide  basis  is  considered  essential  to  further 
improvement  of  long-term  predictions. 

As  the  correlation  between  ionospheric  parameters  and  circuit 
performance  becomes  better  established,  and  as  the  ability  to  predict 
ionospheric  characteristics  in  greater  detail  becomes  available  (e.g. , 
through  the  use  of  topside  sounding),  it  may  prove  useful  to  have  prediction 
methods  at  several  levels  of  sophistication,  the  choice  depending  upon  the 
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knowledge  of  the  circuit  parameters  involved,  the  importance  of  the 
communication  link  and  the  funds  available  to  make  the  predictions. 

16.2.  Short  Term  Predictions 

For  established  circuits  with  an  operating  history,  it  is  often 
advantageous  to  combine  operational  experience  with  predictions  of  the 
type  described  in  this  report  to  optimize  hour-to-hour  frequency  scheduling, 
e.g.  ,  if  frequencies  begin  to  fail  earlier  than  normal  or  if  the  available 
signal-to -noise  ratio  is  lower  than  normal,  the  prediction  may  be  used  to 
anticipate  if  this  trend  is  likely  to  continue.  It  may  also  prove  useful  to 
monitor  selected  frequencies  on  nearby  circuits  to  establish  the  time 
these  frequencies  become  available,  the  available  signal-to-noise  ratio, 
etc.  A  comparison  between  these  observations  and  the  corresponding 
predictions  may  yield  data  to  adjust  the  long  term  predictions  for  the 
circuit  of  interest.  This  approach  is  especially  recommended  where 
operational  history  of  the  circuit  of  interest  is  not  available.  For  high 
priority  operations,  the  use  of  auxiliary  propagation  path  monitors  such 
as  oblique  ionospheric  sounders  can  be  used  to  estimate  the  useful  fre¬ 
quency  limits  and  the  trends  established  from  the  long  term  predictions. 

In  such  applications,  it  may  prove  useful  to  have  predictions  made  for 
both  the  sounder  circuit  and  the  operational  circuit  to  assist  in  the  trans¬ 
lation  of  the  sounder  information  to  expected  trends  on  the  operating 
circuit,  e.  g.  ,  the  useful  frequency  range  on  the  sounder  may  differ  from 
that  of  the  operational  circuit  due  to  differences  in  transmitter  power, 
receiver  sensitivity,  antenna  systems,  etc.  Predictions  for  both  the 
sounder  path  and  the  operational  path  are  expected  to  be  especially  desir¬ 
able  if  it  is  necessary  to  translate  oblique  ionosonde  data  to  operating 
circuits  with  paths  differing  from  that  of  the  sounder. 

For  particularly  difficult  communication  paths,  it  might  prove 
useful  to  prepare  predictions  for  various  solar  activity  levels  and 
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attempt  an  estimate  of  the  appropriate  solar  activity  level  on  a  given 
day.  It  may  be  possible  to  further  refine  the  long  range  predictions  by 
assigning  different  levels  of  effective  solar  activity  to  the  different 
ionospheric  regions,  e.g.  ,  E  and  F2  to  produce  a  family  of  predictions 
from  which  to  make  selections  on  given  days. 

17,  Recommendations 

The  prediction  method  outlined  in  this  report  should  be  extensively 
compared  with  the  observed  performance  cf  high  frequency  sky-wave 
communication  circuits  to  establish  areas  where  modification  or  extension 
of  the  prediction  method  is  required.  It  is  anticipated  that  the  following 

o 

areas  may  merit  particular  attention. 

1.  Nighttime  propagation  losses.  Limited  observations  indicate 
seasonal  and  solar  activity  variations  in  the  nighttime  fields.  The  ability 
to  predict  these  variations  may  be  closely  associated  with  predictions  of 
ionization  levels  in  the  lower  ionospheric  regions,  e.g.  ,  nighttime  D  and 
E  region  ionization  including  sporadic  E. 

2.  Atmospheric  noise  levels  for  directive  antennas.  Current 
atmospheric  noise  levels  are  based  on  observations  made  with  vertically 
polarized  and  essentially  omnidirectional  antennas.  Adjustments  for 
directive  antennas  may  require  a  translation  of  current  noise  maps  into 
noise  source  maps  and  the  use  of  propagation  predictions  to  anticipate 
atmospheric  noise  levels  when  directive  antennas  are  used.  This 
translation  to  noise  sources  could  yield  a  badly  needed  solar  activity 
dependence  of  the  atmospheric  noise  levels  which  could  be  especially 
important  at  the  higher  frequencies. 

3.  Correlation  coefficients  between  the  several  variables  used  in 
the  prediction  method.  Correlation  coefficients  between  signal  levels 
and  noise  levels  are  especially  needed,  but  other  correlation  coefficients 
could  also  materially  add  to  the  prediction  method,  e.g.  ,  sporadic  E  and  F2 
critical  frequencies,  F2  critical  frequencies  and  noise,  etc. 
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4.  Short  term  correlation  coefficients  of  system  loss  and  likelihood 
of  propagation  as  a  function  of  frequency  separation.  A  typical  problem 
is  the  expected  overall  performance  of  a  given  complement  of  frequencies 
or  the  relative  value  of  various  frequencies  in  a  given  complement. 

5.  Practical  antenna  patterns  at  very  low  vertical  angles.  Terrain 
and  atmosphere  may  be  important  factors  in  the  space  wave  radiation  of 
antennas  at  very  low  angles.  Low  angle  paths  are  often  the  only  possible 
paths.  Practical  estimates  of  antenna  gain  at  these  low  angles  could  im¬ 
prove  predictions. 

6.  Ground  reflection  losses.  Although  ground  reflection  losses 
are  not  normally  a  major  factor  in  system  loss  calculations,  the  problem 
of  multiple  hop  propagation  when  the  ground  reflection  involves  extreme 
terrain  conditions  such  as  tropical  rain  forest,  deep  arctic  ice,  or 
mountains,  could  merit  further  investigation.  Side  scatter  from  the 
earth’s  surface  is  part  of  this  problem. 

7.  Ionospheric  scatter  losses.  Ionospheric  scatter  may  be 
important  on  certain  paths  when  regular  ionospheric  reflection  no  longer 
takes  place.  The  combination  of  ionospheric  scatter  and  sporadic-E 
predictions  to  establish  a  floor  in  system  loss  predictions  is  a  possible 
approach. 

8.  Ground  wave  and  troposcatter  losses.  Loss  estimates  for  use 
in  calculations  of  transmitter  and  receiver  site  separation  and  multipath 
between  ground  and  sky  waves  would  be  a  useful  addition  to  computer 
routine. 

9.  Improved  circuit  performance  criteria.  Performance  is  now 
estimated  as  a  function  of  available  sign.al-to-noise  ratio.  Although 
available  signal -to-noise  is  a  basic  criteria,  other  factors  such  as 
multipath  and  doppier  frequency  spreads  may  be  very  significant  especially 
for  high  data  rate  systems.  Comparisons  between  multipath  predictions 
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and  circuit  performance,  and  the  development  of  required  signal-tc-noise 
ratios  both  as  a  function  of  multipath  and  type  of  interfering  noise  could 
be  significant. 

10.  Daily  predictions.  The  prediction  method  of  this  report  is 
based  on  monthly  median  ionospheric  predictions  and  system  losses  and 
the  expected  variation  about  these  values  within  the  month.  The  month 
is  an  arbitrary  time  increment,  and  the  possibility  of  predictions  contin¬ 
uous  in  time  could  merit  investigation. 

The  development  of  D,  E  and  F2-region  predictions  on  this 
ba,sis  could  be  a  first  step.  A  parallel  effort  could  be  the  development 
ofsemiempirica'l  loss  predictions  from  daily  observed  values  in  lieu  of 
the  monthly  median  observed  values. 

11.  Path  geometry.  All  basic  inputs  to  estimating  path  geometry 
require  improvement,  e.g.  ,  F2  critical  frequency.  F2  heigats,  regular 
E  criticals,  sporadic  E,  FI  critical  frequency  and  heights,  plus  other 
appropriate  information  to  describe  the  height  profile  of  ionization.  Ray 
paths  based  on  average  ionospheric  characteristics  could  be  compared 
with  more  sophisticated  ray  tracing  methods  where  ionospheric  charac¬ 
teristics  are  known  in  detail,  e.g.  ,  along  the  path  of  the  topside  sounder, 
to  determine  what  adjustments  would  be  profitable. 

12.  Extension  to  lower  frequencies.  Current  methods  have  a  lower 
limit'  of  3  MHz.  An  extension  to  lower  frequencies  including  a  computer 
program  for  the  medium  frequency  band  is  desirable. 

13.  Inclusion  of  path  antenna  gain  (Gp)  in  lieu  of  the  obvious 
approximation  (Gt  +  Gp)  for  Gp.  [Norton  1959]. 
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Appendix  A.  Theoretical  Antenna  Power  Gain  Equations 


1.  Horizontal  Rhombic  Antenna 


Combining  (8.21  ),  (8.22)  and  (8.23)  and  assuming  the  antenna 
sufficiently  high  for  R  =  600  ohms,  the  gain, g,  of  a  horizontal  rhombic 
is  given  by 


g 


/  ^  TT  \  1  k't  . 

„  „  f  cos  a  sin  **2~  ‘  Ul  *  .  sin  *  Us)'|2 

=  3-n - urns - ■/ 


(A-l) 


•fx  sin2  0  sin2  A  +  X„T  (cos  0  -  cos  A  •  sin  a)2  }• 
L  v  H  J 


where  t  is  the  leg  length  of  the  rhombic  in  meters 
a  is  half  the  obtuse  angle  of  the  rhombic 
0  is  the  azimuthal  angle  of  arrival  measured  counter¬ 
clockwise  from  major  axis  of  the  rhombic 
Uj  =  1  -  cos  A  •  sin  (a  +  0) 

Ua  =  1  -  cos  A  •  sin  (a  -  0). 


2.  Horizontal  Dipole  Antenna 

Using  the  same  approach  as  followed  for  the  horizontal  rhombic 
antenna  the  gain,  g  ,  of  a  horizontal  dipole  of  arbitrary  length,  l,  is 
given  by 


f  .  .  ..  k l 
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where 


R  =  30  |(1  -  cot2  ~)  Cin  (2  k£)  +  4  cot2  ~  Cin  (k£)  (A-3) 

+  2  .  cot  [sA(2kt)  -  2  Si  (k-t)  J  j-  - 

where 
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00 
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3.  Ground  Based  Vertical  Antenna 


For  the  vertical  antenna  we  will  define  P  in  terms  of  the  electric 

r 


field,  E0 , 


and  the  magnetic  field,  H 


0 


by 


P 
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where  E0>  as  given  by  Laitinen  [19  57  J,  is 
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where 


Krb  =  |K*  +  1  +  K2  cos  Uv  -  2b'  ]} 


(A-7) 


R=  30  5  cos  LkfJ  Cin  [  2k-t ]  +  (1  +  cos  [ k-6] )  CinCk-t,] 
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where  t  is  the  height  of  the  vertical  antenna  in  meters. 


4.  Inverted  L  Antenna 
For  the  inverted  L  antenna: 
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.  sin2  Iktl  .  (1  -  u2)} 
u  =  cos  A  sin  0 

v  =  sin  A 

h  =  length  of  vertical  portion  of  inverted  L 

t  =  length  of  horizontal  portion  of  inverted  L 
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R  =  60  {  In  Ck£]  -  Ci  [  kt]  -  —  -  0.4Z3  j 

2k 6 

+  30  [-  0.5  cos  [kh]  Cin  [2kh]  +(1+  cos  [  kh]  )  Cin  Ckh] 

+  sin  [kh]  (0.  5  Si  [2kh]  -  Si  [  kh]  )  } 

0  =  0  is  broadside  to  the  horizontal  portion  of  the  inverted  L. 

The  reflection  factors  ttH  and  KY  are  defined  in  (8.  23b)  and  (8.23a) 
respectively.  The  vertically  polarized  field  of  the  horizontal  portion  of 
the  inverted  L,  is  omitted  in  (A-9).  Its  combination  with  the  vertical 
portion  of  the  inverted  L  to  obtain  the  component  Ng  of  the  radiation 
vector  is  unweildy.  This  is  not  a  serious  omission  since  its  coritribution 
is  small  near  0  =  0. 


5.  Sloping  V  Antenna 


For  sloping  V  antenna: 


g  =  * os{ cos2  a'  (p-- Za  -  05^57 CY‘*W‘-  z=+w=-v=> 

-  y3  (WX.  z4  + wa  .  VA)J J  i-  I - — — - uT 

(y4  (Wa  •  Z3  -  Wi  •  Vs)  -  Ys  {W3  •  Z4  -  Wi  •  V*))  ]2  J  +  [— 27  ’  Zl 


Ui  •  U2 


UCia  •  Zs  K 


Ux  •  Ua 


Us 


—  (w3  (UC45  *  Zs  -  UCss  •  Z4)  +  W4  (UC46  •  Vs 


(w3  (  -  UC^s  «  Vs  +  UCse  *  V4)  +  W4  (UC45  •  Zs  -  UCse  •  Z4))  ]*}  ,  (A-10) 

where 


Ui  =  1  -  cos  Yj 


U3  =  1  -  cos  Ys 
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u2 

=  1  -  COS  f2 

u4 

= 

1  -  cos  f4 

=  Ux  ♦  cos  Ys 

UC36 

= 

U3  •  cos  Ys 

uq* 

-  U2  •  cos  f7 

UC45 

— 

U4  •  cos  fs 

cos  Yx  ■=  sin  A 

sin  A’  + 

cos 

A  cos 

A' 

cos  (0  -  a) 

cos  fa  =  sin  A 

sin  A'  + 

cos 

A  cos 

A* 

cos  (0  +  a) 

cos  f 3  -  -  sin  A 

sin  A’  + 

cos 

A  cos 

A* 

cos  (0  -  a) 

cos  f4  =  -  sin  A 

sin  A*  + 

cos 

A  cos 

A* 

cos  (0  +  ci) 

cos  f 5  =  cos  A 

sin  A*  + 

sin 

A  CGS 

A* 

cos  (0  -  a) 

cos  f 6  =  cos  A 

sin  A'  + 

sin 

COL 

A* 

cos  (0  +  a) 

cos  f7  =  -cos  A 

sin  A'  + 

sin 

A  cos 

A' 

cos  (0  -  a) 

cos  f0  =  -cos  A 

sin  A'  + 

sin 

A  cos 

A’ 

cos  (0  +  a) 

Zx 

=  cos  (k-tUx)  -  1 

z3 

= 

cos  (k-tU3)  -  1 

Za 

=  cos  (ktU2)  -  1 

Z4 

= 

cos  (k-tU4)  -  1 

Vi 

=  sin  (k£Ui) 

v3 

= 

sin  (k-tU3) 

v2 

=  sin  (k-MJ2) 

V4 

= 

sin  (k-tU4) 

Yi 

=  Ux  •  sin  {0  +  a) 

y3 

= 

U3  -  sin  (0  +  a) 

Ya 

=  U2  •  sin  (0  -  a) 

y4 

= 

U4  •  sin  (0  -  a) 

Wi 

=  cos  (Yh  -  2kH  sin  A) 

W3 

= 

cos  (fT  -  2kH  sin  A) 

wa 

=  sin  (Yh  -  2kR  sin  A) 

W* 

= 

sin  (Yy  -  2kH  sin  A) 

H  =  height  of  feed  end  of  sloping  wire  above  ground. 

HT  =  height  of  terminated  end  of  sloping  wire  above  ground 


valid  for  negative  or  positive  slope.  A' 

projection  of  y ,  the  half  angle  betweeu 
the  wires,  onto  the  horizontal  plane, 


A'  =  sin 


.  -1  r  HT  -  H  - 
_  “C  _ 


.  -1  r sin  y  i 

a  =  sin  - rr  i 

Lcos  A'  J 
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where  a  =  [E  +  (m  -  1)  F]  k  sin  A 

E  =  height  of  1st  element  above  ground 

F  =  vertical  separation  between  horizontal  dipoles 
» 

D  =  horizontal  separation  between  feed  points  of  horizontal  dipoles 
N  =  the  number  of  bays  in  the  array 
M  =  the  number  of  elements  in  each  bay 
C  =  distance  to  reflecting  screen 
e  =  0  for  N  even 
1  for  N  odd 

.R  =  rx  -5-  ra  +  .  .  .  +  rM  •  n 
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ri  =  ^ii  **■  &iy  +  •  *  •  ■^■11  +  •  •  •  Rj  m  *  N 

Rjj=  mutual  resistance  of  the  elements  in  the  array  and  is  given  by 

Rj  j  -  -15  cos  [kh]  (  -  2  Ci  A  -  2  Ci  A'  +  Ci  B  +  Ci  B'  +  Ci  C  +  Ci  C'  ) 
+  15  sin  [khj  (2  Si  A  -  Si  A'  -  Si  B  +  Si  B'  -  Si  C  +  Si  C' )  ohms 

A  =  k  { (d2  +  h2)^  +  h  } 

A'  =  k  {  (dc  +  h2)2  -  h] 

B  =  k  { (d2  +  (h  -  if  )^  +  (h  -  l)  } 

B'  =  k  {(d2  +  (h  -  iff  -  (h  -  l)  } 

C  =  k  [  (d2  +  (h  -  if  +  { h  +  l)  } 

C'  =  k  £  (d2  +  (h  +  -t)2)^  -  (h  +  l)  J  , 


where  l  =  length  of  antenna  element 

h  =  horizontal  distance  between  the  ith  and  jth  elements 
d  =  vertical  distance  between  the  ith  and  jth  elements 
R  is  given  by  (A-3). 

7.  Interlaced  Horizontal  Rhombic s 

The  gain  equation  for  two  interlaced  rhombic  antennas  such  that 
they  are  fed  end -fire  is  given  by  (A-l)  where 

Hr  =  fl  +  cos  X  -  Kt  /cos  Yt  (1  +  cos  Z)  -  sin  YT  sin  z\  1 

H  L  H^H  H-’-J 


+  sin  X  -  Kt  {-sinYT  (1  +  cos  Z)  -  cos  YT  sin  Zj  J  , 


where 


X  =  k  [s  -  d  {si 


sin  A  sin  A'  +  cos  A  cos  A'  cos 


*}] 
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Yy  =  Yy  -  2  kH  sin  A 

M  H 

Z  =  X  -  2  kh  sin  A 

S  =  the  horizontal  displacement  between  feed  points  of 
the  two  rhombic s  in  meters 

h  =  the  vertical  displacement  of  the  two  rhombics  in  meters 
H  =  the  height  of  the  lower  rhombic  in  meters 
d  -  (S2  +  h8)*  . 

Also  replace  the  constant,  3.  2,  by  0.  8. 

8.  Sloping  Rhombic  Antenna 
For  the  sloping  rhombic  antenna: 

g  =  •  05  {cos2  A'  [  [Xx  Wx  -  Yi  K„  (cos  Zx  W2  -  sin  Zx  V2)  f 

+  (Xi  V i  -  Yi  Kh  (sin  Zx  W2  +  cos  Zx  V2)  f  ] 

+  {Xa  Wa  -  Ya  Kr  (ccs  Z2  W2  -  sin  Z2  V2)  f 
+  {Xs  Vi  -  Ya  KT  (sin  Z2  W2  +  cos  Z2  V2)  f  j  , 

where 

v  _  sin  [0  -f  a]  sin  [0  -  a] 

Xj  -  -  - 

Ui  u2 
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= 


sin  [0  +  a]  sin  L0  -  ct] 


U3 


ui  =  1  -  sin  A  sin  A'  -  cos  A  cos  A'  cos  (0  4-  a) 


u2  =  1  -  sin  A  sin  A'  -  cos  A  cos  A'  cos  (0  -  a) 


U3  =14  sin  A  sin  A’  -  cos  A  cos  A'  cos  (0  -  a) 


u4  =14-  sin  A  sin  A1  -  cos  A  cos  A' cos  {0  4-  a) 


^2 


cos  A1  sin  A  cos(0  4-a)  -  sin  A*  cos  A  cos  A1  sin  A  cos  (0  -  a] 


ui 


u2 


Y2 


cos  A1  sinA  cos(04-a)4-  sinA*  cos  A 


cos  A'  sin  A  cos  (0  -  a) 


u3 


Zi  =  YH  -  2  JcH  sin  A 

2  T 


Vix  =  1  4-  cos  [k-f.  (uj  4-  u2 )  ]  -  cos  [k-tui]  -  cos  [k-tu2] 
2  34  3  4 


Vi  =  -  sin  [k l  (Ui  4-  u2)  ]  4-  sin  [k^Uj]  +  sin  [k-tu2] 

S  3  4  3  4 


.  -1  rHT  -  HI 

sm  L— — J 


0  and  a  are  the  same  as  defined  for  the  sloping  V. 


-  sin  A!cos  A 


+  sin  A'cos  A 
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Appendix  B,  Tables  of  the  Distribution  of  F2(3000)  MUF 

MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Simspot  Number:  High 

Season:  Winter 

No.  Hemis:  Jan.  ,  Feb.  ,  Nov.  ,  Dec. 
So.  Hemis:  May,  June,  July,  Aug. 

Geog. 

Lat. 

(N  or  S) 

00 

04 

LOCAL  TIME 

08  12 

16 

20 

Upper  Decile 

30° 

1.36 

1. 27 

1.  41 

1.42 

1.  40 

1.  43 

70° 

1.31 

1. 25 

1.  34 

1.30 

1.  16 

1.  34 

60° 

1.26 

1.  23 

1.  24 

1.18 

I.  11 

1.  26 

50° 

1.19 

1.  19 

1.  16 

1.11 

1. 09 

1. 20 

40° 

1.15 

1.  14 

1.  13 

1.09 

1.  09 

1.  14 

30° 

1.22 

1.26 

1.  12 

1.09 

1.  11 

1.  13 

20° 

1.32 

l.  35 

1.  12 

1.  12 

1.  14 

1.  20 

10° 

1.  18 

1.  25 

1.  14 

1.  13 

1.  15 

1. 20 

Lower  Decile 

0 

O 

00 

.62 

.  70 

.  74 

.67 

.  64 

.  73 

o 

© 

.69 

.  74 

.  77 

.  72 

.  72 

.  78 

60° 

.77 

.  78 

.  81 

.80 

.  79 

.  82 

50° 

.83 

.  80 

.  84 

.87 

.  84 

.  86 

o 

O 

.  86 

.  81 

.  87 

.90 

.  87 

.  87 

O 

O 

CO 

.83 

.  76 

.89 

.90 

.  88 

.  86 

20° 

.78 

.  70 

.89 

.89 

.89 

.  83 

O 

O 

H 

.83 

.  76 

.  89 

.90 

.89 

.  84 

Table  B.  1 
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MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Mont  lily  Median) 


Sunspot  Number: 
Season: 

High 

E  quinox 

September 

March,  April,  October 

Geog. 

Lat. 

LOCAL  TIME 

(N  or  S) 

00 

04 

08 

12 

16 

20 

Upper  Decile 

o 

O 

00 

1.46 

1.  37 

1.  35 

1.40 

1. 38 

1. 46 

-o 

o 

o 

1.42 

1.  31 

1.  30 

1.31 

1.  33 

1.  37 

60° 

1.30 

1.  25 

1.  27 

1.24 

1.  25 

1.  24 

<ji 

© 

o 

1. 18 

1.  20 

1.  25 

1.20 

1.  16 

1.  17 

o 

O 

1. 15 

1.  16 

1.  17 

1.  16 

1.  12 

1.  14 

30° 

1.25 

1.  18 

1.  10 

1.  10 

1.  11 

1.  15 

20° 

1.31 

i.  32 

1.  11 

1.  11 

1.  12 

1.  20 

10° 

1.21 

1.  23 

1.  09 

1.20 

1.  14 

1.  23 

Lower  Decile 

80° 

.  66 

.  67 

.  75 

.  66 

.  70 

.  72 

-j 

o 

o 

.67 

.  71 

.  73 

.  70 

.  70 

.  72 

60° 

.69 

.  75 

.  71 

.  71 

.  71 

.  72 

o 

O 

m 

.  73 

.  78 

.  70 

.  72 

.  74 

.  73 

o 

O 

.79 

.  82 

.  75 

.  78 

.  80 

.  82 

OO 

o 

o 

.81 

.  82 

.  87 

.87 

.  87 

.  86 

20° 

.  81 

.  77 

.89 

.92 

.90 

.  85 

10° 

.80 

.79 

.  86 

.90 

.90 

.  82 

Table  B.  2 
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MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number:  High 

Season:  Summer 

No.  Hemis:  May,  June,  July,  Aug. 
So.  Hemis:  Jan.  ,  Feb.  ,  Nov.  ,  Dec 

Geog. 

Lat. 

(N  or  S) 

00 

04 

LOCAL  TIME 

08  12 

16 

20 

Upper  Decile 

80° 

1.  30 

1. 27 

1.  17 

1.15 

1.  23 

1. 24 

70° 

1.22 

1. 22 

1.  20 

1.  18 

1.  21 

1. 23 

60° 

1.  16 

1.  18 

1. 26 

1.21 

1.  19 

1.21 

50° 

1. 14 

1.  15 

1.  30 

1.26 

1.  19 

1.  18 

o 

O 

1.  .14 

1.  14 

1.  30 

1.27 

1.  19 

1.  16 

30° 

1. 16 

1.  15 

1.  25 

1.20 

1.  17 

1.  15 

20° 

1.21 

1.  22 

1.  18 

1.  15 

1.  18 

1.  19 

10° 

1.25 

1.  21 

1.  13 

1.  17 

1.  22 

1.23 

Lower  Decile 


o 

O 

oo 

.73 

.  74 

.  82 

.83 

.79 

.  75 

70® 

.75 

.  75 

77 

.  80 

.  80 

.  77 

60° 

.77 

.  76 

.  74 

.  77 

.  80 

.80 

Ul 

o 

0 

.79 

.  76 

.  73 

.  75 

.  80 

.84 

o 

o 

.80 

.  76 

.  75 

.  75 

.79 

.84 

30° 

.81 

.  76 

.  82 

.81 

.  79 

.  83 

20° 

.81 

.  77 

.  85 

.86 

.  81 

.80 

10° 

.80 

.79 

.86 

.89 

.  85 

.  78 

Table  B.  3 
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MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOU&S 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number: 
Season: 

Medium 

Winter 

No.  Herois:  Jan 
So.  Hemis:  Ma> 

.  ,  Feb.  ,  Nov.  ,  Dec 
r,  June,  July, , Aug. 

Geog. 

Lat. 

(N  or  S) 

00 

04 

LOCAL  TIME 

08  12 

16 

20 

Upper  Decile 

80° 

1.45 

1. 39 

1. 44 

1.40 

1. 33 

1:  45 

9 

o 

t"- 

1.39 

1.  31 

I.  37 

1.32 

1.  29 

1.41 

o 

o 

vO 

1.33 

1.  24 

1. 25 

1.21 

1.  22 

1.  33 

U1 

o 

9 

1.30 

1.  19 

1.  14 

1.  15 

1.  16 

1.  29 

9 

o 

1.27 

1.  17 

1.  12 

1. 14 

1.  14 

1.  28 

o 

o 

1.30 

1.  31 

1.  16 

1. 18 

1.  18 

1.  32 

20° 

1.33 

1.  38 

1.  17 

1.22 

1.  26 

1.  40 

10° 

1.21 

1.26 

1.  14 

1.13 

1.  15 

1.  23 

Lower  Decile 

80° 

.76 

.  78 

.  68 

.67 

.  62 

.  70 

70° 

.79 

.81 

.  74 

.70 

.  73 

.  73 

60° 

.82 

.  83 

.79 

.75 

.  80 

.  76 

o 

O 

m 

.84 

.  82 

.  83 

.81 

.  84 

.  78 

9 

o 

.83 

.  81 

.  85 

.86 

.  86 

.79 

o 

o 

.  73 

,  76 

.  85 

.85 

.  85 

.  78 

20® 

.  74 

.  71 

.  85 

.83 

'.82  ' 

.  76 

O 

O 

H 

.  77 

.  69 

.  87 

.86. 

.  85  ' 

.  78 

Table  B.  4 


MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number: 
Season: 

Medium 

Equinox 

March,  April, 

September, 

October 

Geog. 

Lat. 

(N  or  S) 

00 

04 

LOCAL ‘ 

08 

riME 

12 

16 

20 

Upper  Decile 

80° 

1.45 

1.31 

1.  27 

1.28 

1.  30 

1.47 

o 

O 

1.41 

1.  22 

1.  23 

1.26 

1.  26 

1.  38 

60° 

1.35 

1.  17 

1. 20 

1.23 

1.  18 

1. 29 

50° 

1.28 

1.  15 

1.  17 

1.21 

1.  13 

1.  20 

40° 

1.22 

1.  16 

1.  16 

1. 18 

1.  12 

1.  17 

3(  0 

1.22 

1.22 

1.  15 

1. 17 

1.  14 

1.23 

o 

O 

C4 

1.32 

1.  30 

1.  13 

1.15 

1.  17 

1.  37 

10° 

1. 18 

1.39 

1.  11 

1. 13 

1.  20 

1.23 

Lower  Decile 

00 

o 

o 

.64 

.61 

.  73 

.  74 

.  74 

.  67 

70° 

.68 

.  71 

.  77 

.  74 

.  78 

.  70 

60° 

.70 

.  75 

.  80 

.72 

.  78 

.  73 

o 

e 

.73 

.  77 

.  81 

.  74 

.  76 

.  75 

0 

O 

.75 

.  78 

.  82 

.VS 

.  76 

.  76 

30° 

.77 

.  76 

.  82 

.83 

.  78 

.  72 

20° 

.75 

.  73 

.  84 

.87 

.  81 

.69 

10° 

.79 

.  68 

.  86 

.89 

.  84 

.  80 

Table  B.  5 
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MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number: 

Medium 

No.  Hemis:  May, 

June, 

July, 

Aug. 

Season: 

Summer 

So.  Hemis:  Jan.  , 

Feb. 

,  Nov. 

,  Dec. 

Geog. 

Lat. 

LOCAL  TIME 

(N  or  S) 

00 

04 

08 

12 

16 

20 

Upper  Decile 

00 

o 

o 

1.27 

1. 23 

1. 20 

1. 18 

1.  25 

1.  23 

70° 

1.23 

1.  19 

1.  19 

1.  17 

1.  17 

1.  19 

60° 

1.20 

1.  18 

1.  19 

1.  17 

1.  14 

1.  17 

50° 

1. 17 

1.  19 

1.  21 

1.  17 

1.  15 

1.  16 

40° 

1. 17 

1.22 

1.  23 

1.  13 

1.  17 

1.  17 

30° 

1.20 

1.  30 

1.  22 

1. 19 

1.  19 

1.  18 

20° 

J.26 

1.  38 

1.  17 

1.23 

1.  23 

1.  28 

10* 

1.26 

1.44 

1.  11 

1.28 

1. 28 

1.  22 

Lower  Decile 

o 

O 

00 

.82 

.  80 

.  82 

.85 

.  80 

.  79 

70° 

.83 

.82 

.  79 

.82 

.  82 

.  82 

60° 

.83 

.  82 

.  77 

.79 

.82 

.  83 

50° 

.81 

.  81 

.  76 

.77 

.  81 

.  82 

40° 

.78 

.  78 

.  75 

.78 

78 

.  78 

30° 

.  77 

.  73 

.  75 

.79 

.  77 

.  74 

20° 

.77 

.  69 

.  78 

.82 

.  78 

.  73 

10° 

.79 

.  63 

.  84 

.85 

.  81 

.  77 

Table  B.  6 
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MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number: 
Season: 

Lov; 

Winter 

No.  Hemis:  Jan. 
So.  Hemis:  May, 

Feb.  , 
June, 

Nov.  , 
July, 

Dec. 

Aug. 

Geog. 

Lat. 

(N  or  S) 

00 

04 

LOCAL  TIME 

08  12 

16 

20 

Upper  Decile 

o 

O 

00 

1.44 

1.  34 

1. 45 

1.32 

1. 33 

1.40 

o 

o 

t- 

1.37 

1. 29 

1.  38 

1.23 

1. 24 

1. 35 

60° 

1.30 

1.  24 

1.  27 

1.  15 

1.  17 

1. 30 

50° 

1.25 

1.  21 

1.  16 

1.  12 

1.  12 

1. 25 

o 

o 

1.23 

1.  20 

1.  13 

1.11 

1.  11 

1. 23 

30° 

1. 28 

1.  30 

1.  15 

1.  17 

1.  15 

1.28 

20° 

1.34 

1.  37 

1.  19 

1.20 

1.  24 

1.32 

10° 

1.27 

1.  38 

1.  18 

1.  15 

1.  14 

1.  20 

Lower  Decile 

80° 

.60 

.  65 

.69 

.  72 

.  68 

.  67 

o 

O 

.68 

.  71 

.  75 

.76 

.  75 

.  70 

60® 

.74 

.  76 

.  80 

.80 

.  82 

.  73 

50° 

.79 

.  78 

.  83 

.85 

.84 

.  76 

o 

O 

.81 

.79 

.  85 

.87 

.89 

.  77 

o 

o 

.81 

.  74 

.  86 

.32 

.  85 

.  78 

o 

O 

C\] 

.78 

.  67 

.  87 

.  75 

.  77 

.79 

10° 

.71 

.  70 

.  88 

.86 

.87 

.79 
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MUF  EXCEEDED  0.  10  AND  0.  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number: 
Season: 

— .  .  k  - 

Low 

Equinox 

March, 

April, 

September , 

October 

Geog. 

Lat. 

(N  or  S) 

00 

04 

LOCAL  TIME 

08  12 

16 

20 

Upper  Decile 

0 0 
o 

o 

1.42 

1.  32 

1. 29 

1.26 

1.  33 

1. 48 

70° 

1.38 

1.  25 

1. 25 

1.23 

1.  26 

1.  40 

60° 

1.32 

1. 21 

1. 22 

1.20 

1.  20 

1.  31 

50° 

1.26 

1.  19 

1. 20 

1.  18 

1.  16 

1.  26 

9 

o 

1.22 

1.  20 

1.  19 

1.  16 

1.  16 

1. 25 

30° 

1.22 

1.  26 

1.  18 

1.15 

1.  16 

1.  28 

20° 

1.30 

1.  32 

1.  16 

1.  14 

1.  18 

1.  33 

O 

O 

H 

1.23 

1.40 

1.  13 

1.  13 

1.  19 

1.  16 

Lower  Decile 

80° 

.67 

.  72 

.  74 

.  73 

.  80 

.  65 

70* 

.  70 

.  75 

.  76 

.  74 

.  82 

.  69 

60* 

.73 

.  78 

.  80 

.  75 

.  81 

.  73 

50° 

.75 

.  80 

.  81 

.76 

.  81 

.  76 

40° 

.  77 

.  81 

.  81 

.  77 

.  80 

.  78 

30° 

.78 

.  80 

.  82 

.  78 

.  81 

.  74 

20° 

.  77 

.  75 

.  83 

.81 

.  83 

.69 

10° 

.  76 

.  66 

.  86 

.89 

.  86 

.  75 
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MUF  EXCEEDED  0.  10  AND  0,  90  OF  THE  HOURS 
(Ratio  of  MUF  to  Monthly  Median) 


Sunspot  Number: 

Low 

No.  Hemis:  May, 

June , 

July, 

Aug. 

Season: 

Summer 

So.  Hemis  :  Jan.  , 

Feb.  , 

Nov. 

,  Dec. 

Geog. 

Lat. 

LOCAL  TIME 

(N  or  S) 

00 

04 

08 

12 

16 

20 

U  pper  Decile 

00 

o 

o 

1.26 

1.  24 

1.  15 

1.  17 

1. 21 

1.  22 

70° 

1. 22 

1.  18 

1.  14 

1.15 

1.  16 

1.  18 

60° 

1.18 

1.  17 

1.  14 

1.15 

1.  14 

1.  15 

50° 

1.  17 

1.  20 

1.  15 

1. 16 

1.  14 

1.  15 

o 

O 

1.  17 

1.  25 

1.  17 

1.17 

1.  15 

1.  16 

30° 

1.  18 

1. 30 

1.  17 

1.20 

1.  19 

1. 20 

to 

o 

o 

1. 20 

1. 34 

1.  14 

1.24 

1.  22 

1.  23 

10° 

1.20 

1. 37 

1.  12 

1.30 

1.  27 

1.  20 

Lower  Decile 

o 

O 

00 

.68 

.79 

.  84 

.87 

.  85 

.  76 

o 

O 

r- 

.70 

.  81 

.  83 

.  86 

.  86 

.  77 

60° 

.7? 

.  84 

.  83 

.84 

.  86 

.  81 

50° 

.75 

.  85 

.  82 

.83 

.  85 

.  84 

o 

O 

■<* 

.79 

.  85 

.  80 

.82 

.  83 

.  85 

30° 

.79 

.  82 

.  78 

.80 

.  81 

.  80 

20° 

.77 

.  78 

.  77 

.79 

.  79 

.  73 

10° 

.74 

.  75 

.  80 

.83 

.  82 

.  69 
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Appendix  C.  Tables  of  the  Distribution  of  Transmission  Loss 

EXPECTED  EXCESS  SYSTEM  LOSS  ABOV  E  QUASI-MINIMUM  (dB) 
(Whiter  -  Paths  <2500  km) 


01-04  LMT 

04-07  LMT 

G.  M. 
Lat. 

Med. 

st 

% 

Med. 

st 

Su 

0-40 c 

9.0 

4.  0 

9.0 

9.0 

4.0 

7.  6 

40-45° 

9.4 

4.3 

9.0 

9.3 

4.3 

8.  3 

45-50° 

9.9 

4.7 

9.  1 

9.6 

4.6 

9.0 

50-55° 

10.3 

5.  1 

9.2 

10.  0 

5.0 

9.7 

55-60° 

11.0 

5.3 

10.  0 

13.  5 

6.  7 

9.6 

60-65° 

16.  C 

8.  0 

13.  5 

25.0 

12.  7 

13.  0 

65-70° 

15.2 

7.7 

14.  6 

24.6 

13.  5 

13.  2 

70-75° 

.12.4 

6.3 

9.4 

17.  3 

8.9 

15.  2 

75-80° 

11.0 

5.6 

9.4 

15.  6 

7.  7 

8.8 

07-10  LMT 

10-13  LMT 

0-40° 

9.0 

4.0 

7.  6 

9.0 

4.  0 

6.4 

40-45° 

10.  1 

4.6 

8.6 

9. 1 

4.  5 

7.  1 

45-50° 

li.  2 

5.2 

9.6 

9.2 

5.  1 

7.8 

50-55° 

12.2 

5.9 

10.  7 

9.3 

5.7 

8.7 

55-60° 

15.6 

8.2 

14.  6 

10.4 

5.0 

10.  6 

60-65° 

23.0 

12.3 

23.7 

12.6 

6.8 

20.  5 

65-70° 

21.8 

11.8 

22.  5 

11.3 

6.  0 

22.  0 

70-75° 

17.  6 

9.9 

14.3 

9.9 

5.4 

13.9 

75-80° 

15.2 

8.4 

10.2 

10.  5 

6.3 

10.  7 
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T 


EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI-MINIMUM  (dB) 


(Winter 

-  Paths 

<  2500  km) 

13  -  16  LMT 

16  -  19  LMT 

G.  M. 
Lat. 

Med. 

St 

Su 

Med. 

St 

Su 

0-40° 

9.  0 

4.  0 

6.  4 

9.  0 

4.  0 

7.  6 

40-45° 

9.  1 

4.  2 

6.  2 

9.  6 

4.  6 

8.  1 

45-50° 

9.  2 

4.  4 

6.  5 

10.  2 

5.  2 

8.  6 

50-55° 

9.  4 

4.  6 

6.  9 

10.  8 

5.  8 

9.  1 

55-60° 

9.  8 

4.  8 

7.  2 

11.  7 

6.  5 

9.  0 

60-65° 

11.  0 

5.  8 

8.  7 

15.  8 

8.  3 

14.  1 

65-70° 

10.  5 

5.4 

8.  2 

13.  6 

7.  8 

11.  3 

70-75° 

9.  5 

4.  8 

7.  5 

11.  5 

6.  5 

10.  5 

75-80° 

9.  5 

4.  7 

6.  7 

10.  4 

5.4 

8.  6 

19  -  22 

LMT 

22  -  01  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

9.  0 

40-45° 

9.  7 

4.  6 

7.9 

9.3 

4.  2 

9.  1 

45-50° 

10.  5 

5.  2 

8.  3 

9.  6 

4.  5 

9.  2 

50-5$ 

11.  8 

5.  8 

8.  7 

9.  9 

4.  8 

9.  3 

55-6<f 

12.  0 

5.  4 

10.  6 

11.  0 

5.  1 

9.  5 

60-65° 

18.  8 

8.  9 

17.  2 

15.  6 

7.  2 

9.  9 

65-7  Cf 

16.  0 

7.  8 

18.  6 

14.  0 

6.  8 

11.  5 

7  0-7$ 

12.  0 

5.  8 

14.  8 

12.  0 

6.  0 

9.  0 

75-8(5* 

10.  5 

4.  9 

11.  5 

10.  0 

5.  4 

8.  5 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI-MINIMUM  (dB) 
(Winter  -  Paths  >  2500  km) 


01-04  LMT 

04 

-07  LMT 

G.  M. 
Lat. 

Med. 

st 

Su 

Med. 

s;. 

Su 

0-40° 

9.  0 

4.  0 

9.  0 

9.  0 

4.  0 

7.  6 

40-45° 

9.  3 

4.  2 

9.  1 

9.  3 

4.  3 

7.  8 

45-50° 

9.  6 

4.  4 

9.  3 

9.  6 

4.  6 

8.  0 

50-55° 

9.  9 

4.  6 

9.  5 

10.  0 

5.  0 

8.  2 

55-60° 

10.  6 

4.  7 

9.  6 

11.  4 

5.  5 

8.  3 

60-65° 

14.  3 

5.  7 

11.  4 

16.  0 

6.  8 

9.  5 

65-70° 

14.  5 

6.  5 

10.  2 

14.  0 

6.  3 

10.  9 

70-75° 

11.  5 

3.  4 

9.  7 

13.  5 

6.  5 

7.  8 

75-80° 

10.  5 

5.  1 

9.  2 

11.  7 

5.  9 

8.  1 

07-10  LMT 

10 

-13  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

6.  4 

40-45° 

9.  4 

4.  2 

9.  0 

9.  0 

4.  2 

7.  3 

45-50° 

9.  8 

4.  5 

10.  4 

9.  1 

4.  4 

8.  2 

50-55° 

10.  3 

4.  8 

11.  9 

9.  1 

4.  6 

9.  2 

55-60° 

11.  7 

5.  6 

12.  4 

9.  6 

5.  2 

10.  4 

60-65° 

14.  9 

7.  3 

14.  1 

11.  0 

5.  2 

15.  2 

65-70° 

16.  0 

8.  3 

14.  2 

10.  9 

4.  4 

15.  8 

70-75° 

13.  3 

6.  8 

11.  2 

9.  8 

4.  8 

11.  2 

75-80° 

11.  5 

5.  7 

10.  ? 

9.  8 

4.  8 

9.  2 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI- MINIMUM  (dB) 

(Winter  -  Paths  >  2500  km) 


13-16  LMT 

16-19  LMT 

G.  M. 
Lat. 

Med. 

Su 

Med. 

st 

su 

0-40° 

9.  0 

4.  0 

6.  4 

9.  0 

4.  0 

7.  6 

40-45° 

9.  0 

4.  0 

6.  5 

9.  3 

4.  3 

7.  7 

45-50° 

9.  1 

4.  1 

6.  7 

9.  6 

4.  7 

7.  8 

50-55° 

9.  1 

4.  1 

6.  9 

10.  0 

5.  1 

7.  9 

55-60° 

9.  4 

4.  6 

7.  4 

10.  7 

5.  5 

7,  6 

60-65 

9.  7 

4.  1 

8.  3 

13.  0 

6.  7 

8.  0 

65-70° 

9.  5 

5.  1 

7.  7 

12.  0 

6.  2 

7.  3 

70-75° 

9.  3 

4.  7 

7.  1 

10.  6 

5.  4 

7.  4 

75-80° 

8.  8 

4.  4 

7.  0 

9.  9 

5.  0 

7.  5 

19-22  LMT 

22-01  LMT 

0-40 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

9.  0 

40-45° 

9.  5 

4.  3 

7.  6 

9.  4 

4.  3 

9.  2 

45-50° 

10.  0 

4.  6 

7.  6 

9.  5 

4.  6 

9.  4 

50-55° 

10.  5 

4.  9 

9.  6 

10.  4 

4.  9 

9.  6 

55-60° 

10.  8 

5.  0 

7.  9 

12.  5 

5.  3 

9.  8 

60-65° 

13.  5 

5.  6 

8.  9 

19.  0 

7.  8 

13.  8 

65-70° 

13.  5 

6.  1 

8.  5 

18.  6 

8.  6 

11.  1 

70-75° 

11.  0 

5.  1 

8.  6 

13.  8 

6.  5 

9.  2 

75-80° 

9.  8 

4.  8 

8.  2 

10.  6 

5.  0 

9.7 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI-MINIMUM  (d3) 


(Equinox  - 

Paths  < 

2500  km) 

01-04  LMT 

04- 

07  LMT 

G.  M. 
L-at. 

Med. 

Sa 

Med. 

s-t. 

Su 

0-40° 

9.  0 

4.  0 

9.  0 

9.  0 

4.  0 

9.  0 

40-45° 

9.  5 

4.  5 

10.  0 

9.  6 

4.  4 

11.  5 

45-50° 

10.  1 

5.  0 

11.  1 

10.  3 

4.  8 

14.  1 

50-55° 

10.  7 

5.  6 

12.  2 

11.  0 

5.  2 

16.  6 

55-60° 

11.  4 

5.  7 

17.  6 

13.  4 

6.  4 

22.  0 

60-65° 

15.  7 

7.  7 

30.  3 

20.  2 

9.  5 

29.  3 

65-70° 

15.  5 

8.  1 

28.  0 

21.  0 

11.  1 

31.  0 

70-75° 

17.  3 

7.  0 

21.  7 

20.  0 

13.  8 

20.  8 

75-80° 

10.  4 

6.  1 

15.  5 

11.  5 

7.  5 

18.  7 

07-10  LMT 

10-13  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

6.  4 

40-45° 

10.  6 

5.  3 

9.  8 

10.  0 

4.  7 

9.  0 

4^-50° 

12.  3 

6.  6 

12.  0 

11.  0 

5.  4 

11.  6 

50-55° 

14.  0 

8.  0 

14.  3 

12.  0 

6.  2 

14.  2 

55-60° 

a  6.  5 

8.  3 

15.  3 

14.  0 

7.  6 

18.  3 

60-65° 

26.  0 

14.  0 

23.  4 

18.  0 

10.  6 

33.  0 

65-70° 

30.  4 

18.  2 

26.  9 

17.  5 

10.  0 

27.  9 

70-75° 

20.  6 

12.  8 

20.  2 

14.  1 

8.  8 

18.  9 

75-80° 

16.  4 

9.  7 

14.  4 

12.  8 

7.  5 

13.  6 

Table  C.  5 


EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI- MINIMUM* dB} 
(Equinox  -  P?ths  <  2500  km) 


13 

-16  LMT 

16- 

19  LMT 

G.  M. 

Lat. 

Med. 

% 

Su 

Med 

% 

Su 

0-40° 

9.  0 

4.  0 

6.  4 

9.  0 

4.  0 

7,  6 

40-45° 

9.7 

4.  5 

8.  9 

10.  3 

5.  0 

11.  3 

45-50° 

10.  4 

5.  0 

11.  4 

11.  6 

6.  0 

15.  0 

50-55° 

11.  2 

5.  6 

13.  9 

13.  0 

7.  0 

18.  7 

55-60° 

11.  6 

5.  6 

15.  5 

13.  8 

7.  5 

20.  2 

60-65° 

16,  2 

8.  3 

19.  2 

18.  0 

10.  3 

27.  0 

65-70° 

13.  8 

7.  0 

18.  0 

15.  0 

8.  4 

24.  0 

70-75° 

11.  6 

6.  2 

14.  2 

13.  0 

7.  2 

18.  0 

75-80° 

10.  0 

5.  4 

12.  0 

11.  4 

6.  2 

14.  1 

19-22  LMT 

22-01  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

9.  0 

40-45° 

10.  0 

4.  8 

10.  0 

10.  3 

4.  7 

10.  0 

45-50° 

11.  0 

5.  6 

12.  5 

11.  6 

5.  4 

11.  0 

50-55° 

12.  0 

6.  4 

15.  0 

13.  0 

6.  2 

12.  0 

55-60° 

14.  5 

7.  7 

19.  5 

15.  1 

7.  4 

13.  3 

60-65° 

19.  9 

11.  3 

29.  0 

24.  0 

13.  0 

26.  7 

65-70° 

19.  0 

11.  3 

28.  8 

22.  7 

11.  2 

17.  5 

70-75° 

15.  0 

8.  6 

22.  0 

16.  0 

8.  0 

16.  5 

75-80° 

11.  4 

6.  4 

20.  6 

12.  3 

6.  3 

15.  7 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI- MINIMUM  (dB) 
(Equinox  -  Paths  >  2500  km) 


01-04  LMT 

04-07  LMT 

G,  M 
Lat. 

Med. 

Su 

Med. 

Su 

0-40° 

9.  0 

4  0 

9.  0 

9.  0 

4.  0 

7.  6 

40-45° 

9.  3 

4.  1 

10.  0 

9.  3 

4.  1 

8.  5 

45-50° 

9.  6 

4.  2 

11.  0 

9.  6 

4.  2 

9.4 

50-55° 

9.9 

4.  4 

12.  1 

9.  9 

4.  3 

10.  3 

55-60° 

10.  4 

4.  5 

13.  2 

11.  0 

4.  6 

10.  6 

60-65° 

12.  9 

5.  7 

15.  5 

14.  2 

5.  9 

10.  8 

65-70° 

12.  7 

5.  7 

14.  3 

14.  6 

6.  6 

10.  6 

70-75° 

10.  8 

4.  9 

13.  1 

11.  9 

5.  3 

9.  8 

75-80° 

10.  0 

4.  8 

11.  0 

10.  2 

4.  6 

9.  0 

07-10  LMT 

10-13  LMT 

o 

i 

o 

0 

9.  0 

4.  0 

7.  6 

vO 

o 

4.  0 

6.4 

40-45° 

9.  6 

4.  2 

8.  3 

9.  6 

4.  4 

7.9 

45-50° 

10.  2 

4.  5 

9.  0 

10.  3 

4.  9 

9.4 

50-55° 

10.  8 

4.  8 

9.  7 

11.  0 

5.  4 

11.  0 

55-60° 

12.  2 

5.  7 

9.  8 

12.  8 

6.  3 

11.  2 

60-65 

16.  6 

7.  9 

11.  4 

18.  0 

9.  4 

15.  7 

65-70° 

16.  7 

7.  7 

13.  8 

17.  1 

9.  0 

13.  4 

70-75  ° 

13.  1 

6.  1 

10.  9 

14.  7 

7.  7 

12.  4 

75-80° 

11.  7 

5.  7 

10.  6 

12.  0 

6.  5 

12.  2 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI- MINIMUM  (dB) 


( Equinox 

-  Paths  > 

2500  km) 

13-16  LMT 

16- 

19  LMT 

G.  M. 
Lat. 

Med. 

St 

Su 

Med 

St 

0-40° 

9.  0 

4.  0 

6.  4 

9.  0 

4.  0 

7.  6 

40-45° 

9.  4 

4.  3 

8.  8 

10.  0 

4.  5 

9.  0 

45-50° 

9.  8 

4.  6 

11.  2 

11.  0 

4.  9 

10.  4 

50-55° 

10.  3 

5.  0 

13.  7 

12.  1 

5.  4 

11.  9 

55-60° 

10.  8 

5.  1 

14.  2 

13.  0 

5.  7 

12.  0 

60-65° 

14.  0 

6.  5 

17.  2 

16.  5 

7.  5 

14.  2 

65-70° 

12.  5 

5.  7 

17.  3 

14.  0 

6.  5 

13.  3 

70-75° 

10.  5 

4.  9 

15.  9 

12.  2 

5.  8 

12.  5 

75-80° 

9.  6 

5.  0 

14.  8 

11.  0 

5.  5 

11.  0 

19-22  LMT 

22-01  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

9.  0 

40-45° 

9.9 

4.  3 

8.9 

9.  6 

4.  2 

10.  6 

45-50° 

10.  8 

4.  6 

10.  2 

10.  3 

4.  4 

12.  3 

50-55° 

11.  7 

5.  0 

11.  6 

11.  0 

4.  7 

14.  0 

55-60° 

14.  4 

6.  2 

12.  4 

12.  4 

5.  1 

17.  5 

60-65° 

19.  0 

8.  5 

16.  7 

16.  9 

7.  5 

22.  1 

65-70° 

18.  5 

8.  5 

16.  0 

16.  3 

7.  1 

17.  5 

70-75° 

14.  6 

6.  9 

13.  8 

12.  6 

5.  6 

16.  3 

75-80° 

11.  3 

5.  5 

13.  7 

10.  6 

4.  9 

15.  6 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI- MINIMUM  (dB) 
(Summer  -  Paths  <  2500  km) 


01 

-04  LMT 

04-07  LMT 

G.  M. 
Lat. 

Med. 

St 

Su 

Med. 

St 

Su 

0-40° 

9.  0 

4.  0 

9.  0 

9.  0 

4.  0 

7.  6 

40-45° 

9.  8 

4.  4 

9.  8 

9.  6 

4.  3 

9.  9 

45-50° 

10.  6 

4.  8 

10.  6 

10.  2 

4.  7 

12.  2 

50-55° 

11.  5 

5.  3 

11.  4 

10.  8 

5.  1 

14.  6 

55-60° 

12.  2 

5.  5 

17.  8 

11.  7 

5.  5 

16.  1 

60-65° 

15.  5 

7.  5 

24.  5 

13.  8 

5.  8 

22.  7 

65-70° 

13.  9 

6.  5 

22.  1 

13.  2 

6.  2 

21.  8 

70-75° 

11.  4 

5.  4 

15.  6 

12.  0 

5.  8 

15.  0 

75-80° 

11.  2 

5.  7 

12.  8 

11.  7 

6.  0 

11.  3 

07-10  LMT 

10-13  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

6.  4 

40-45° 

9.  8 

4.7 

9.  1 

9.7 

4.  5 

7.  2 

45-50° 

10.  6 

5.  5 

10.  7 

10.  5 

5.  0 

8.  1 

50-55° 

11.  4 

6.  3 

12.  3 

11.  3 

5.  5 

9.  0 

55-60° 

13.  0 

7.  2 

15.  6 

12.  0 

6.  0 

10.  6 

60-65° 

18.  0 

10.  7 

26.  0 

15.  0 

7.  8 

18.  8 

65-70° 

15.  2 

9.  2 

26.  7 

13.  5 

6.  7 

19.  5 

70-75° 

11.  8 

6.  8 

18.  2 

12.  0 

6.  2 

12.  0 

75-80° 

10.  2 

5.  3 

16.  8 

11.  3 

5.  7 

9.  5 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI-MINIMUM  (dB) 
(Summer  -  Paths  <  2500  km) 


13- 

■16  LMT 

16-19  LMT 

Ct.  m. 
Lat. 

Med. 

% 

Su 

Med. 

St 

Su 

0-40° 

9.  0 

4.  0 

6.  4 

9.  0 

4.  0 

7.  6 

40-45° 

9.  8 

4.  5 

8.  1 

10.  6 

5.  1 

8.  6 

45-50° 

1.0.  7 

5.  0 

9.  8 

12.  3 

6.  2 

9.  7 

50-55° 

11.  6 

5.  5 

11.  6 

14.  0 

7.  3 

10.  8 

55-60° 

12.  0 

5.  8 

13.  5 

15.  5 

8.  5 

11.  0 

60-65° 

13.  8 

6.  2 

19.  9 

18.  1 

9.  1 

14.  2 

65-70° 

12.  8 

5.  7 

16.  6 

14.  4 

7.  2 

13.  8 

70-75° 

11.7 

5.  5 

13.  3 

13.  2 

6.  9 

11.  1 

75-80° 

11.  2 

5.  3 

12.  1 

12.  8 

7.  0 

10.  2 

19-22  LMT 

22-01  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

9.  0 

40-45  ° 

10.  5 

4.9 

9.  5 

10.  3 

4.7 

10.  0 

45-50° 

12.  1 

5.9 

11.4 

11.  6 

5.  4 

11.  0 

50-55° 

13.  6 

6.  9 

13.  4 

13.  0 

6.  1 

12.  0 

55-60° 

15.7 

7.  8 

14.  1 

14.  0 

6.4 

15.  0 

60-65° 

18.  5 

10.  1 

20.  5 

17.  9 

8.9 

21.  1 

65-70° 

16.  7 

7.  9 

21.  9 

16.  2 

7.  3 

19.  0 

70-75° 

13.  8 

7.  6 

17.  8 

13.  5 

6.  8 

16.  0 

75-80° 

13.  2 

7,  2 

11.  8 

12.  8 

6.  8 

12.  7 

Table  C.  10 


159 


- - ^ 


EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI-MINIMUM  (dB) 
(Summer  -  Paths  >  2500  km) 


01- 

04  LMT 

04- 

07  LMT 

G.  M. 
Lat. 

Med. 

sc 

Su 

Med. 

St 

Su 

0-40° 

9.  0 

4.  0 

9.  0 

9.  0 

4.  0 

7.  6 

40-45° 

9.  7 

4.  4 

9.  1 

9.  8 

4.  4 

9.  1 

45-50° 

10.  4 

4.  8 

9.  2 

10.  6 

4.  9 

10.  6 

50-55° 

11.  2 

5.  2 

9.4 

11.  4 

5.  4 

12.  2 

53-60° 

11.  8 

5.  4 

9.  6 

12.  8 

6.  2 

13.  0 

60-65° 

14.  5 

6.  7 

9.  8 

17.  5 

8.  8 

16.  8 

65-70° 

13.  4 

6.  1 

10.  0 

15.  3 

7.  4 

16.  7 

70-75° 

10.  9 

4.  8 

8.  9 

13.  6 

6.  5 

11.  9 

75-80° 

10.  8 

5.  3 

8.  2 

12.  8 

6.  0 

10.  0 

07-10  LMT 

10-13  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

6.  4 

40-45° 

9.  7 

4.  5 

8.  1 

9.  3 

4.  2 

6.  9 

45-50° 

10.  4 

5.  0 

8.  6 

9.  6 

4.  4 

7.  5 

50-55° 

11.  2 

5.  6 

9.  2 

10.  0 

4.7 

8.  1 

55-606 

12.  7 

6.  5 

9.7 

10.  4 

4.  9 

9.  2 

60-65° 

16.  6 

9.  3 

13.  8 

11.  6 

5.  6 

13.  1 

65-70° 

14.  9 

8.  2 

14.  5 

10.  9 

5.  4 

13.  1 

70-75° 

11.  9 

5.  9 

11.  1 

10.  3 

5.  0 

10.  1 

75-80° 

11.  2 

5.  5 

10.  1 

10.  1 

4.  8 

8.  6 
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EXPECTED  EXCESS  SYSTEM  LOSS  ABOVE  QUASI-MINIMUM  (dB) 
(Summer  -  Paths  >  2500  km) 


13-16  LMT 

16- 

19  LMT 

G.  M. 
Lat. 

Med. 

St 

Su 

Med. 

St 

Su 

0-40° 

9.  0 

4.  0 

6.  4 

9.  0 

4.  0 

7.  6 

40-45° 

9.  7 

4.  5 

7.  6 

9.  4 

4.  2 

8.  1 

45-50° 

10.  4 

5.  0 

8.  8 

9.  8 

4.  5 

8.  6 

50-55° 

11.  2 

5.  6 

10.  1 

10.  2 

4.  8 

9.  1 

55-60° 

11.  7 

5.  9 

12.  3 

10.  6 

5.  1 

9.7 

60-65° 

13.  2 

6.  8 

17.  4 

12.  2 

6.  1 

11.  3 

65-70° 

12.  3 

6.  3 

14.  0 

10.  7 

5.  4 

10.  4 

70-75° 

11.  2 

5.  5 

12.  2 

10.  3 

5.  1 

9.  0 

75-80° 

11.  1 

5.  4 

11.  4 

10.  2 

5.  0 

8.  2 

19-22  LMT 

22-01  LMT 

0-40° 

9.  0 

4.  0 

7.  6 

9.  0 

4.  0 

9.  0 

40-45° 

10.  4 

4.  9 

8.  2 

9.  8 

4.  4 

9.  6 

45-50° 

11.  9 

5.  8 

8.  9 

10.  6 

4.  9 

10.  3 

50-55° 

13.  4 

6.  7 

9.  6 

11.  5 

5.  4 

11.  0 

55-60° 

13.  5 

6.  8 

11.  8 

12.  2 

5.  7 

13.  4 

60-65° 

14.  5 

7.  1 

17.  2 

14.  7 

7.  2 

18.  4 

65-70° 

13.  4 

6.  4 

15.  2 

13.  5 

6.  4 

15.  1 

70-75° 

11.  8 

i 

5.  8 

12.  3 

11.  7 

5.  3 

12.  1 

75-80° 

11.  3 

5.  6 

9.9 

11.  4 

5.  5 

10.  3 
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Appendix  D:  Tables  of  Prediction  Errors  in  Excess 

System  Loss  Above  Quasi-Minimum  (dB) 

Winter  Months 


Geomag. 
Lat  • 

01-07  LMT 

°.u 

°p 

07-13  LMT 

°.u 

G>1 

0-40° 

5.60 

.  60 

1.  63 

3.  10 

.80 

1.89 

40-45° 

5.61 

.  62 

1. 65 

3.  13 

.  62 

2.  16 

45-50° 

5.71 

.  65 

1.66 

3.  28 

.  68 

2.25 

50-55° 

5.80 

.  68 

1.  67 

3.  30 

.  78 

2.  30 

55-60° 

5.75 

.  88 

1.88 

3.  58 

1.26 

4.  50 

60-65° 

6.75 

2.48 

2.05 

5.  15 

2.48 

5,48 

65-70° 

6.82 

2.  50 

2.  16 

4.  93 

2.49 

5.31 

70-75° 

5.93 

1.  85 

2.92 

4.  10 

1.80 

2.96 

75-80° 

5.83 

2.  08 

2.08 

3.  63 

1.27 

2.  07 

13-19  LMT 

19-01  LMT 

0-40° 

3.  10 

.  60 

1.68 

5.  60 

.60 

1.66 

40-45° 

3.20 

.85 

1.  78 

5.  62 

.63 

1.  67 

45-50° 

3.32 

.87 

1.80 

5.  64 

.  65 

1.  70 

50-55° 

3.  56 

1.00 

2.  30 

5.  65 

.72 

1.  75 

55-60° 

4.08 

1.25 

2.  90 

5.  66 

1.00 

1.80 

60-65° 

7.43 

2.  65 

3.62 

5.97 

1.35 

3.  65 

65-70° 

5.03 

1.  73 

3.  66 

5.91 

1.09 

4.  10 

70-75° 

4.  12 

1.25 

6.  06 

5.  80 

.67 

2.80 

75-80° 

3.45 

1.23 

2.05 

5.  63 

.63 

2.  75 
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Prediction  Errors  in  Excess  System  Loss 
Above  Quasi-Minimum  (dB) 


Summer  Months 


Geomag. 

Lat. 

01-07  LMT 

^5 

07-13  LMT 

°p 

a>l 

0-40° 

5.60 

.  60 

1.  65 

3.  10 

.60 

1.  61 

40-45° 

5.63 

.  65 

1.  68 

3.  21 

.62 

1 .  66 

45-50° 

5.  70 

.  78 

1.81 

3.30 

.  73 

1.  80 

U1 

o 

i 

oi 

0 

5.86 

1. 07 

2.  57 

3.  47 

1. 00 

2.  06 

55-60° 

5.  88 

1.  35 

3.  01 

3.  94 

1.  52 

2.  67 

60-65° 

6.  56 

2.  57 

4.84 

5.  71 

2.86 

4.  73 

65-70° 

6.42 

2.  05 

4.  51 

4.  87 

2.  16 

6.  38 

70-75° 

5.88 

1.40 

2.76 

3.  67 

1.26 

2.95 

75-80° 

5.85 

1.25 

2.01 

3.  55 

1.24 

2.90 

13-19  LMT 

19-01  LMT 

0-40° 

3.  10 

.  60 

1.68 

5.  60 

.  60 

1.66 

40-45° 

3.25 

.73 

1. 77 

5.  63 

1.25 

1.  85 

45-50° 

4.  52 

1.22 

1.92 

6.  10 

2.28 

2.  34 

50-55° 

4.  58 

1.  71 

2.30 

6.  15 

3.29 

3.05 

55-60° 

5.  15 

2.  27 

2.62 

8.  02 

3.31 

3.40 

60-65° 

8.62 

3.78 

4.  36 

8.  00 

3.85 

3.40 

65-70° 

5.22 

2.09 

3.22 

7.  80 

2.25 

2.95 

70-75° 

4.40 

1.72 

2.  61 

6.  56 

1.99 

2.83 

75-80° 

4.  25 

1. 65 

2.  57 

6.  32 

1.72 

2.  05 
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Prediction  Errors  in  Excess  System  Loss 
Above  Quasi-Minimum  (dB) 


Equinox  Months 


Geomag. 

Lat. 

ap 

01-07  LMT 

a.u 

a‘l 

ap 

07-13  LMT 

°.u 

0-40° 

5.  60 

.  60 

1.  63 

3.  10 

.  60 

1.  60 

40-45° 

5.68 

.  62 

1.71 

3.  28 

.  61 

1.  76 

45-50° 

5.69 

.  65 

2.02 

3.34 

1.  00 

2.  42 

50-55° 

5.92 

1.24 

2.  61 

3.  75 

1.  52 

3.  08 

55-60° 

6.38 

1.  82 

4.  09 

4.  40 

1.81 

3.  28 

60-65° 

9.05 

3.28 

6.45 

5.  70 

7.  65 

7.35 

65-70° 

9.92 

4.  54 

7.  26 

7.  75 

4.  85 

8.  62 

70-75° 

7.  10 

3.  05 

4.98 

4.  93 

2.49 

5.48 

75-80' 

5.87 

1.35 

3.  62 

4.  35 

2.09 

3.08 

13-19  LMT 

19-01  LMT 

0-40° 

3.  10 

.  60 

1.  68 

5.  60 

.  60 

1.66 

40-45° 

3.  32 

.  72 

1.  77 

5.  63 

.63 

1.68 

45-50° 

3.87 

1.25 

2.  38 

5.  70 

.72 

1.95 

50-55° 

4.98 

1.  73 

2.  84 

5.  78 

.91 

2.24 

55-60° 

5.  63 

1.  75 

2.  85 

6.27 

1.35 

2.96 

60-65° 

7.85 

3.  75 

5.  06 

8.30 

2.28 

5.30 

65-70° 

6.45 

2.69 

4.  50 

7.81 

2.29 

5.71 

70-75° 

4.  85 

1.  73 

2.93 

6.  32 

1.  35 

3.45 

75-80° 

3.89 

1.  25 

2.  61 

5.  92 

1.  07 

2.  78 
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APPENDIX  E:  TABLES  OF  THE  PREDICTION  ERRORS  IN  THE 

ATMOSPHERIC  RADIO  NOISE  LEVELS 
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Numerical  Coefficients  for  use  in  Evaluation  of  the 
Standard  Deviation  of  the  Lower  Decile  Value  of  the  Atmospheric  Radio  Noise  (cr0L ) 
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